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I. INTRODUCTION 



A. BACKGROUND 

The switching function as performed in microwave circuitry utilized in radar 
applications is essential for the performance of the circuit objective. In the majority of 
practical cases, switching must be done by semiconductor devices, since mechanical 
switches are ruled out by requirements of speed and allowable system weight and size. 
Therefore the semiconductor switch is a vital component in most microwave circuit 
applications. 

Despite its practical importance, relatively little research attention has been given 
to fundamental aspects of the switching problem, in comparison with, for example, the 
problem of microwave amplifier or filter design. Therefore it is the purpose of this 
thesis to examine some of the basic features of the switching problem. A general 
requirement in the design of a switching circuit is the conversion of the set of given 
switch performance specifications into a set of device requirements for the necessary 
semiconductors to be used in the switch. It is known that a device figure of merit can 
be defined which provides a unifying factor in this conversion process. Some progress 
in the definition and use of switch figures of merit has been made by earlier workers, 
but few practical procedures have been developed, and further work is needed in this 
area. It is the purpose of this thesis to investigate the characterization of several 
aspects of switch performance as a function of device figure of merit, as well as to 
investigate the important aspect of the selection and use of devices so as to maximize 
the power-handling capacity of the switch. 

The fundamental specifications governing microwave switch selection and/or 
design are: 

1. Low insertion loss in on-state 

2. High isolation in off-state 

3. Broad bandwidth 

4. Sufficient power capability 

These characteristics are interactive, a change in one affecting the others. The 
microwave switches considered here are those which employ semiconductor devices 
(p-i-n diodes and FET's) operated in two alternate bias states. The device is 
characterized by a pair of two-terminal impedance values Zj and Z2 in the respective 
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bias states. The parasitic capacitance associated with semiconductor switching devices 
necessitates some form of impedance compensation in order to achieve acceptable 
insertion losses for the switch. Consequently the theoretical analysis of switching 
circuitry has been oriented toward the study of two-state impedances embedded in 
lossless impedance transformers [Refs. 1.2,3]. These researches on switching theory 
have introduced specific invariant parameters related to the device impedance in its two 
states, under lossless impedance transformation. 

Kurokawa and Schlosser [Ref. 1], showed that diodes incorporated in a one-port 
reflective switch can be characterized by a quality factor parameter, Q~ : 

q 2 = — 1 — ^ (eqn 1.1) 

R 1 R 2 

where Rj is the resistive component ofZj, and 
R 2 is the resistive component of Z 2 * 

The numerical value of Q is unchanged when the impedances are subjected to a lossless 
impedance transformation. 

Atwater and Sudbury [Ref. 2] showed that this Q had the character of a figure of 
merit, serving as an indicator of switching performance for semiconductors embedded 
in variety of circuit configurations, and not restricted to the original reflection-type 
circuit of Kurokawa and Schlosser. They followed the research of Kurokawa and 
Schlosser and assumed that Zj and Z 2 had been transformed to real (resistive) values 
in both states. Equation 1.1 becomes: 

q 2 = t R r . R ?l ( e qn L2 ) 

R 1 R 2 

From equation 1.2, we have the relationships: 

R> is a function of Q and R->. and 
R-> is a function of Q and R,. 

Knowing the value of Q and given that Rj, R 2 are known, the switching performance 
is fully specified. 
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Kawakami [Ref. 3] found a different form of switching invariant parameter. He 
defined this parameter as a semiconductor figure of merit, M: 



|Z r Z 2 | 



(eqn 1.3) 



The. parameter M assumes values between 0 and 1, and is also unaltered by lossless 
impedance transformation. Kawakami also showed that when a device was 
incorporated in a symmetrical two-port junction as a switch, its transmission 
coefficients in the on and off states, S 2 ^(l) and S 2 ^( 2 ), respectively, are related by 
equation: 



|S 21 (2)| 



1S 2 i( 1)1(1-M) 
1-M El 1-2 S 21 (l)|] 



(eqn 1.4) 



Equation 1.4 fully characterizes the switching performance of the device, with 
parameter M as a figure of merit. It may inverted to yield the form: 



M 



iS 2I (2)i-|S 21 (l)| 

|S 21 (2)| + |S 21 (1)|[1-2|S 21 (2)|] 



(eqn 1.5) 



This expression is based on the assumption that the switching device in its two states 
has the properties: 



r i 

r 2 



0 (matched; Z = Zq ) 
r max = M 



Equation 1.5 presents a useful tool for the selection of a semiconductor device to 
satisfy low power switching applications. 1 

A fundamental theorem was derived by Hines [Ref. 4], to predict the power that 
can be switched by a single-diode switch when it is connected in a iumped-element 
circuit such as to resonate the parasitic reactance of the switch at a single frequency. 
Hines showed that the power which can be switched by the diode, P sw , is given by: 



J In the absorption switch assumed in equation 1.4, the diode absorbs the incident 
energy in the switch off-state. 
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0-5[V max I m ax^ 



(eqn 1.6) 



P 



sw 



where V max is a diode breakdown voltage, and 

I ma x a maximum forward current rating. 

Hines conjectured that P represents a theoretical limit for power switching by a 
given device in an impedance-matched circuit. 

In many applications, microwave switching circuits typically require multiple 
diodes in cascade connection to achieve adequate levels of off-state isolation. For 
example, radar and phased-array modules need a double-throw switch of high isolation 
to protect the receiver front end. 

B. SCOPE OF THE THESIS RESEARCH 

For the application to microwave switching, the proposed semiconductor figures 
of merit are tested for different types of switching circuitry, to find the best parameter 
to use in a switch selection and/or design problem in order to minimize insertion loss in 
the on-state, and maximize isolation in the off-state. Additional objectives are to 
determine, for the designed switch, what are the power dissipation and heat removal 
problem in the individual semiconductor components. 
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II. CIRCUIT ANALYSIS 



In the present chapter we generalize the circuit representation of microwave 
switches of our interest. Following this, a suitable circuit model may be selected for 
the particular type of switch required. The simplest representation is that of an N- 
port, with scattering matrix description. The number of ports is determined by the 
number of semiconductor devices and the type of switch (single-throw or multiple- 
throw). For the semiconductor devices, there are three possible ways of connection to 
our circuit model, in shunt or in series or a combination of both. Because of heat- 
sinking considerations we select the shunt connection and use a one-port network to 
represent the two-terminal semiconductor device. 

A. CLASSIFICATION OF MICROWAVE SWITCHES 

1. Single Pole Single Throw Switches (SPST) 

The single pole - single throw switch is a switch which has one input port and 
one output port. In this topology we consider our switch as a 'box' and it is a two-port 
box network shown in Figure 2.1(a). For an ideal switch, all the power is transferred to 
the load at the output port in its on-state, and no power transferred to the load in its 
off-state. 

2. Single Pole N Throw Switches (SPNT) 

For these types of switches, there is an input port and N output ports. Their 
switching circuit models are a three-port box network for single pole double throw 
switches, shown in Figure 2.1(b), and a (N+l)-port box network for single pole N- 
throw switches, shown in Figure 2.1(c). In microwave switch application, we are 
interested in single pole double throw switches as mentioned in the introductory 
chapter. For SPDT switches, when one arm of the switch is in the on-state, the other 
arm is in the off-state, and vice versa. In general, we may consider a SPDT switch as 
two SPST switches with their input arms tied together and appiv two different bias 
states to the semiconductor devices in the different arms of the switch. 



2 Single pole is a line-over-ground circuit. 
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B. SCATTERING MATRIX REPRESENTATION OF MICROWAVE 
SWITCHES 



We can use the properties of the S-matrix [Refs. 5,6: pp. 112, 50], defined in 
Figure 2.2, and assume that our switches are ideal switches with matched load at both 
input and output ports. We can represent the scattering matrices for both on and off 
states of our switches as follows. 

1. Single Pole Single Throw Switch 



ON STATE: [S] 



0 1 
1 0 



OFF STATE: [S] 



1 0 
0 1 



2. Single Pole Double Throw Switch 
ARM (2) ON: [S] = 



0 1 0 
1 0 0 
0 0 0 



ARM (3) ON: [S] - 



0 0 1 
0 0 0 
1 0 0 



C. SWITCHING CIRCUIT MORPHOLOGIES 



1. Number of Semiconductor Devices 

In an SPST switch or one arm of SPDT switch there is at least one 
semiconductor device. For improvement of load isolation, the switch may require 
additional semiconductor devices. In general, a SPNT switch may use M 
semiconductor devices. 



M = pN (eqn 2.1) 

where M is the number of semiconductor devices, 

N is the number of throws, and 
p is a positive integer. 
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a. SPST 



3 PORT 2 . 

BOX Z 



b. SPOT 






portl 

'L 



N PORT 
NETWORK 





V 2 port2 



«e-Vn^ 
port n 



b l = s ll a l + s 12 a 2 + + S]_ n a n 

b n = s nl a l + S n2 a 2 + + s nn a n 



Figure 2.2 Scattering Parameters. 
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2. Switch "BOX" Topologies 

We can define the switching circuit as a lossless embedding network with one 
input port, N output ports, and M ports to which the semiconductor devices are 
coupled. The total number of ports of the embedding network is given by: 



Example of switching circuits. 

1. SPST switch with 1 semiconductor device. Figure 2.3(a) 

Number of ports = (1 + 1 + 1) 

= 3 ports 

2. SPST switch with 2 semiconductor devices. Figure 2.3(b) 

Number of ports = (1 + 1 + 2) 

= 4 ports 

3. SPDT switch with 2 semiconductor devices. Figure 2.3(c) 

Number of ports = (2+1 + 2) 

= 5 ports 

D. SWITCHING SEMICONDUCTOR DEVICE CHARACTERIZATION 

The semiconductor device in shunt connection can be represented by a one-port 
network, as stated above. 3 The device itself in the switching circuit may be 
characterized by: 

1. Two-states of impedance [Ref. 2]. Zj ( i = 1,2 ) or 

2. Two-states of reflection coefficient [Ref. 1]. T- ( i = 1,2 ). 

The relations between impedance and reflection coefficient are: 



Number of ports 



(N+l + M) 



(eqn 2.2) 




(eqn 2.3) 




(eqn 2.4) 



3 The device is in shunt to ground in the single pole circuit. 
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1 



2 



i 1 




a. SPST 1 DEVICE 




b. SPST 2 DEVICES 




c. SPDT 2 DEVICES 

Figure 2.3 Example or' Switching Circuits. 
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Figure of merit of the semiconductor device is defined based on the device 
characterization chosen. 

1. Figure of Merit, M 
From equation 1.3 



M ' 



|Z r Z 2 | 

IZi-Zj* 



Kawakami postulated that the fundamental characterization was in terms of 
two specific states of . He assumed the following conditions: 



r j = 0 ( Z is matched; Z = Zq ) 



A lossless matching network is assumed to be used to transform the device 
impedance to the matched condition in this state. In the second bias state of the 
device, the reflection coefficient presented by the transformed device impedance is then: 



|r 2 | = |F max ( ( this is equal to the figure of merit, M, of the devices ) 

However, this type of switch is a reflection switch. 4 

2. Quality Factor, Q 

Kurokawa and Schlosser postulated that the fundamental characterization was 
in terms of two states of impedance and their real parts (resistive). They defined a 
quality factor Q , shown in equation 1.1, as the figure of merit of the semiconductor 
devices. 

From the parameters defined above we can characterize the p-i-n diode and 
specify a diode quality factor or the figure of merit for the device used in microwave 
switching circuits. Diode quality factor, Q^iode’ ^ or a switching diode is equal 

to the switching quality factor, as may be shown: 

From equation 1.1 

, (R, - Ro ) 2 (X, - X ,) 2 

Q*- = — 1 

RjR-. 



4 The T- = 0 state, absorbs energy and is not suitable for high power switching 
circuits. 1 
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For a p-i-n diode, a simple circuit model is assumed in which the two 
impedance states are: 

Zj = Rj forward bias 

Z 2 = R 2 + jX ? reverse bias 

If we assume that Fq = R 2 = R ( for both bias states ) 




1 

( RCco ) 2 

Thus Kurokawa-Schlosser's Q reduces to the conventional diode Q: 



Q 2 = Q 2 diode ( ec l n 2 ’ 6 ) 

3. Comparison of Figure of Merit 

We can find the relationships between the figure of merit, M, and the quality 
factor, Q, as follows: 

Kurokawa-Schlosser, quality factor, from equation 1.1 



0 2 = IZ|- 2/ 

R 1 R 2 

2 _ (R|-R 2 ) 2 + (X r X 2 ) 2 

RjR 2 

Kawakami, figure of merit, from equation 1.3 



|M| 2 



iz 1 -z 2 i 2 
|Z 1 + z 2 "| 2 



|M| 2 



(Ri - R 2 )“ - (Xi - x 2 >- 

(Rj + R 2 ) 2 + (Xj -X 2 > 2 



1 

1 - |M| 2 



, IZi-Z 2 l 2 
IZ, + Z 2 1 2 



(eqn 2.7) 



(eqn 2.8) 
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iz, + z 2 v - |Z, - Z 2 | 2 
|Z| + Z,*| 2 




4|M| 2 




, |M| 2 
4 *-) 

MMp 



(eqn 2.9) 



|M| 2 




(eqn 2.10) 



If Q 2 > > 1 

|M| = \.lj (eqn 2.11) 

E. TRANSMISSION LINE SECTION REPRESENTATION OF SWITCHING 
CIRCUIT MODEL 

From an N-port box representation, the semiconductor devices are in shunt to 
ground and have two alternate impedance values. For practical microwave integrated 
circuit (MIC) applications, it is convenient to use transmission line sections in the 
embedding network of the switches 

Figure 2.4, shows the circuit model of SPST and SPDT switches with one 
semiconductor device in each arm. To optimize the switch characteristics for low 
insertion loss and high isolation, the impedance values of the devices are kept 

unchanged at the particular frequency in each state. The length and characteristic 
impedance of the transmission line section are vaned to achieve the optimum switching 

characteristics in a particular frequency interval. 
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1 



* ? — 2 

Zo I Zo 

J SEMICONDUCTOR 
| DEVICE 



S?ST 1 DEVICE 



1 



* — 3- 




Zo 



w 3 * 

— : * 2 

ZO 

SEMICONDUCTOR 

DEVICE 




SEMICONDUCTOR 

DEVICE 



3 



SPDT 2 DEVICES 

0 = Transmission Line Electrical Length 

Zo = Transmission Line Characteristic Impedance 



Figure 2.4 Transmission Line Model of Switch Networks. 
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F. APPROACHES TO COMPUTATION OF CIRCUIT POWER- 

DISSIPATION 

Using the Touchstone program, we can design a switch to have optimum 
switching characteristics. But the semiconductor devices in the optimized circuit may 
not be able to handle the amount of power dissipated due to heating of the devices. 
For effective switching design, it is advisable to calculate the maximum pow r er 
dissipation in each device before construction of the switch. 

1. ABCD Matrix Representation of a Two-port Network 

From Figure 2.5 the ABCD matrix is defined in equations 2.12 and 2.13 



V j — AV 2 B 1 2 

h = cv 2 + di 2 

Equations 2.12 and 2.13 have the matrix form. 



' v l' 




'A 


B 




' V 2‘ 


.*1. 




_C 


D 




- l 2- 


< N 
> 




‘A 


B' 


-1 


' V f 


h 




_C 


D. 




LiJ 





'D 


-C' 




’V 




_-B 


A. 




Ihi 



(eqn 2.12) 
(eqn 2.13) 

(eqn 2.14) 



(eqn 2.15) 



The utility of the general circuit parameter matrix lies in the fact that cascaded 
two-ports can be represented by the chained product of element matrices [Ref. 5: p 
110 ]. 



V 




A i 


B i' 




A 2 


B 2 _ 






B n" 




' V n+r 


_h . 




_ c i 


D i 




_ C 2 


O 

K) 




. C n 


D n. 




Jn+lJ 



(eqn 2.16) 



Thus the general circuit parameter matrix representation of cascaded two- 

ports is the matrix product of general circuit parameter matrices of the individual two- 
ports. as shown in Figure 2.6. 

2. ABCD Matrix Representation of the Switching Circuit Model 

In the circuit model we have transmission line sections and semiconductor 
devices combined. Lumped-circuits can be characterized by ABCD two-port matrices 
[Ref. 7: p. 18S], as follows: 
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Figure 2.5 ABCD Matrix Model. 





Figure 2.6 N Two-Port in Cascade with Equivalent Matrix. 
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2.1 Series impedance, as shown in Figure 2.7 (a) 



A 


b' 




'1 


z' 


.C 


D. 




.0 


1 . 



(eqn 2.17) 



2.2 Shunt admittance, as shown in Figure 2.7 (b') 



'A 


B' 




'l 


0 ' 


_C 


D. 




Y 


1_ 



(eqn 2.18) 



2.3 Lossless transmission line section, as shown in Figure 2.7 (c) 

"A B 
.C D 

0 is an electrical length of a transmission line section 



Cos0 jZoSinG 
jYASin0 Cos0 



(eqn 2.19) 



2.4 Lossy transmission line section, as shown in Figure 2.7 (d) 



'A B" 




Cosh(y 1) 


ZoSinh(Y 1 ) 


_C D. 




. Y 0 Sinh(Y 1) 


Cosh(Y 1) 



1 = Physical length of a transmission line section, 

y = a + jj3 

a = transmission line attenuation constant. 5 
P = wavenumber. 

3. Calculation of Power in a Cascade Circuit 
From Figures 2.5, and 2.6 [Ref. 8: p. 393], 



(eqn 2.20) 



Z in 


AZ l + B 
CZ L + D 


(eqn 2.21) 


I. 


V 2 


( eqn 2.22) 


1 

N 

(1-4 

1 

N 




v i 


7- 

in v 
Z + Z- 8 


(eqn 2.23) 


V* 


= (Ajlj) * (CjVj) 


(eqn 2.24 ) 



5 a =0 for a lossless transmission line section. 
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Figure 2.7 ABCD Matrix Representation of Devices. 
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Vj +1 = (DjVj) ' (Bjlj) 

j = 1,2,3, n 

Power dissipated in each two-port network. 

P k = 0.5^e[(V k I k *)-(V k+1 I k+1 ’ > )] 

k = 1,2,3, ,n 

Power at load 

P L = 0.5^e[V n+1 I n+1 ^] 

" P in-I p k 
k= i 



(eqn 2.25) 



(eqn 2.26) 



(eqn 2.27) 



(eqn 2.28) 
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III. MICROWAVE SWITCHING DIODE 



The use of the semiconductor diode as a switching element in microwave circuits 
is based on the difference between the diode impedance under reverse and forward bias. 
The diode appears as a very small impedance under forward bias and a very large 
impedance under reverse bias. For a shunt-connected diode as a microwave switch, in 
its on state most of the power is delivered to the load, and the diode looks like a large 
impedance (reverse bias state). On the other hand, in its off state, the power is 
reflected back to the source, or dissipated in the diodes, and the diodes look like a 
small impedance (forward bias state) in shunt to ground. In this chapter we will 
introduce some characteristics of the p-i-n diode, its behavior under forward and 
reverse bias states, the equivalent circuit, and power handling capability. 

A. THE P-I-N DIODE STRUCTURE 

In a p-i-n diode, the semiconductor wafer has a heavily doped P-region and a 
heavily doped N-region separated by a layer of high resistivity material that is nearly 
intrinsic. The thickness of the high resistivity layer lies in the range between 1.5 to 150 
microns (Ref. 9; p. 5]. The electrical contacts are made to the two heavily doped 
regions. 

In practice, there are no materials without any impurities. A practical p-i-n 
diode consists of extremely high resistivity P or N-type material between low resistivity 
(highly doped) P and N-regions. The lightly doped, high resistivity layer is called 7t-type 
for P material, and v-type for N material. The I-region of a p-i-n diode can consist of 
either v or 7t-type material. Figure 3.1(a) shows a P-v-N diode structure. If the I- 
region is of sufficiently high resistivity, with few impurity atoms that will be ionized 
the depletion zone will extend throughout the I-region and include a small 
penetration into both the P and N-regions. Because of heavy doping in the P and N- 
regions the depletion zone will not extend very far into them, and the width of the 
depletion zone will be essentially equal to the I -layer width, W. Figure 3.1(b) shows 
the structure of P-7t-N diode, the depletion zone width is approximately equal to the 
width of the intrinsic region. Most p-i-n diodes use v material [Ref. 7: p. 41] for the 
I -layer. 
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B. THE ELECTRICAL EQUIVALENT CIRCUIT MODEL 

The p-i-n diode used in microwave application has a simplified equivalent circuit 
under forward and reverse bias states as shown in Figure 3.2. 

1. Forward Bias Equivalent Circuit 

Under the forward bias state, both holes and electrons are injected into the I- 
region, and the mobile carriers form a conductive plasma whose resistivity depends on 
the carrier density obtained from the injection process. The forward bias resistance, Rf 
may be expressed as: 

Rf = Rj + R s p (eqn 3.1) 

where R- = resistance due to the I region, and 

R s f = resistance of P and N-regions and their ohmic contacts. 

For an abrupt junction diode, the resistance due to the I-region can be 
expressed as follow: 



R: 






2 if H av 



(eqn 3.2) 



where W = I-region thickness, 

if = DC forward bias current, 
p av = the average carrier mobility, and 

= the carrier lifetime in the intrinsic region. 

In the shunt configuration switching application the forward bias resistance, 
R-, is the major factor of power dissipation in the diode. Since Rj is an inverse function 
of a DC bias current, if , we can increase if to minimize R' . For low voltage p-i-n 
diodes, an if of 10 mA, and for high voltage p-i-n diodes 150 to 200 mA is sufficient 
to make Rf equal to R g f [Ref. 9: p. 8]. 

2. Reverse Bias Equivalent Circuit 

Under the reverse biased state, the I-region is depleted of carriers and the p-i-n 

diode appears as a constant capacitance to a microwave signal. The dissipation losses 
occur in the ohmic contacts and the resistances of the P and N-regions. R f is used to 
represent these losses in the equivalent circuit model. Due to the capacitive reactance 
at microwave frequency, the circuit impedance of the p-i-n diode in its reverse bias 
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a. P-l'-N IMPURITY PROFILE 



I (IT) 




b. P-f-H IMPURITY PROFILE 

Figure 3.1 Profiles for the two P-I-N Diodes Types. 



R . 



J 




FORWARD BIAS REVERSE BIAS 

Figure 3.2 P-I-N Diode Chip Equivalent Circuit. 
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state is much larger than the impedance in the forward bias state. The capacitance, Cj 
of the p-i-n diode is determined by the diode size and I-region thickness as shown in 
equation 3.3 



C0£ r 7t 

4W 



(eqn 3.3) 



where £q = free space permittivity, 

£ f = the relative dielectric constant of the material, 

D = the junction diameter, and 
W = I- region thickness. 

This equation is valid for all reverse bias voltage magnitudes greater than the 
minimum voltage for complete I-layer depletion. Figure 3.3 shows equation 3.3 
graphically for typical available p-i-n diode I-regions. 

3. P-I-N Diode Cutoff Frequency 

Hines [Ref. 4], defined a figure of merit for p-i-n diodes to relate the switching 
effectiveness in terms of switching cutoff frequency, f cs . The Hines figure of merit is 
given by 



‘cs 



27tC-VRfR r 



(eqn 3.4) 



The figure of merit, f cs is one of the parameters used in microwave switching 
design, and for Rp = R r , f cs /f is equal to the conventional diode Q. Knowing Cj, Rp, 
and R we can estimate the operating frequency range of the selected p-i-n diode from 
equation 3.4 

C. POWER HANDLING 

In the off-state of the microwave switch, the RF incident power is reflected back 
to the source or dissipated in the p-i-n diode. This dissipation causes the temperature 
of the junction to rise. For each diode there is a maximum temperature that it can 
withstand. The other factor which limits the amount of power handling is the peak 
voltage that appears across the junction, which can cause damage to the junction. In 
the p-i-n diode specification sheets these important parameters are available for the 
designer to estimate the power handling capability of the diode in the designed 
switches. 
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1. Maximum Voltage across the P-I-N Diode 

The instantaneous voltage across a p-i-n diode junction is the algebraic sum of 
the RF voltage and DC bias voltage. If the junction voltage exceeds the device’s 
breakdown voltage, the p-i-n diode is in the avalanche region, and the p-i-n diode acts 
as a low resistance causing high negative current flow and rapid temperature rise. 
Figure 3.4 shows the condition for a high RF voltage drive with reverse bias. The 
maximum voltage across the junction occurs when the diode is in the reverse bias state. 
The maximum instantaneous voltage applied to the junction is defined as: 

^max ~ ^ ^Fpeak + ^bias ( ec l n 

where V^- as is the DC reverse bias impressed on the diode. 

White [Ref. 7: p. 155], suggested the conservative limit of the safe maximum 
voltage across the junction as: 



V 



max 



0.5 V b 



(eqn 3.6) 



where V b is the p-i-n diode breakdown voltage. A circuit designer may choose a 
different safety factor within an acceptable margin of safety for the designed switch. 

2. Maximum Junction Temperature 

The junction temperature of a p-i-n diode is determined by the ambient 
temperature of the circuit and the power dissipated in the diode. The dissipated power 
may be pulsed or CW. In either case, junction temperature must be kept below the 
maximum allowable value. This maximum temperature is dependent on the diode 
material (semiconductor), and the materials used in mounting and connecting to the 
semiconductor. The recommended maximum junction temperature for safe operation is 
200° C [Ref. 10: p.290]. 

a. CW Power Dissipation 

In the case of a continuous or a long-duration pulse. Tie temperature limit 
,s determined by the initial temperature of the diode materials, the temperature of the 

heat sink, and thermal resistance between the wafer and the heat sink [Ref. 10: p.290]. 
0j c is the total thermal resistance unit in °C /watt. The junction temperature, Tj, for a 
given power dissipated, and the ambient temperature, T a is given by the equation: 



- T a + p d e jc 



(eqn 3.7) 
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JUNCTION CAPAClfANCf C.lj.fl 




Figure 3.3 P-I-N Diode Cj vs. I-Region Diameter and Thickness. 




Figure 3.4 Junction Voltage at High RF and Reverse Bias. 
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The value of 0j c is determined by the physical volume of the junction, 
material and geometry of the package. The maximum value of 0;. for most glass 
packages can be as much as 600°C/watt. A large-area, high capacitance chip mounted 
on copper heat sink can have extremely low 0j c of 5° C/watt or less [Ref. 9: p. 10]. 
b. Pulse Power Dissipation 

For a short duration pulse, the temperature limit is determined by the 
thermal time constant x ^ of the active region. 

t th = OjHC (eqn 3.8) 

where 0j is the junction region thermal resistance, this value is available in the 
specification sheet or can be calculated from the diode physical size. 



6j - 



w 



K t h A 



(eqn 3.9) 



where W = the width of the I-region, 

= the thermal conductivity of the material, 

A = the junction area, and 

HC = the thermal capacity of the active region, also available in the 
specification sheet or can be calculated from: 



HC = C p pWA (eqn 3.10) 

where Cp = the specific heat of the material, and 
p = the density of the material. 

For a pulse width of 0.25 or less, the pulse width is small compared 
with the thermal time constant. All the thermal energy is absorbed in the active region 

area, not spreading to the surrounding. The temperature rise of the junction, AT, in 
units of °C is expressed by: 

P H At 

AT = -S— (eqn 3.11) 

HC 
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where At is the pulse duration. The junction temperature 6 Tj for a short pulse width 
from the above assumption is the sum of the ambient temperature and the temperature 
rise of the junction. 

Tj = T a + AT (eqn 3.12) 

For pulse widths of O.oOt^ to 3r t ^ the diode junction temperature follows 
an exponential curve given by the equation 

Tj = T a + P d 0j [1 - exp(-t/T th ) ] (eqn 3.13) 

Note that equation 3.13 contain 0j, the junction thermal resistance and not 0j c the 
total thermal resistance. For cooling between pulses, the junction temperature follows 
an exponential decay with the same time constant [Ref. 9: p. 10]. 

D. TABLE OF TYPICALLY AVAILABLE P-I-N DIODES 

To illustrate typical p-i-n diodes available for microwave switch design, Table 1 
lists a number of typical diode chips from Microwave Associates, Inc. and their 
approximate electrical parameters. The symbols used in this table are already defined in 
this chapter. 



6 Neglecting the heat flow from the junction. 
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TABLE 1 

TYPICAL PARAMETERS OF P-I-N DIODES 



MODEL 

NUMBER 


MA- 

4P102 


MA- 

4P103 


MA- 

4P202 


MA- 

4P203 


MA- 

4P303 


MA- 

4P404 


MA- 

4P504 


Minimum 

Vb(Volt) 


50 


30 


100 


100 


200 


300 


500 


Maximum 
Cj @ Vr(10) 
(pFXVolt) 


0.5 


.15 


.05 


.15 


.15 


.20 


.20 


Maximum 
Rs @ III 10) 
(ohm)(ma) 


2.5 


1.5 


2.0 


1.5 


1.5 


0.6 


0.6 


Typical 

Carrier 

Lifetime(jiS) 


.100 


.010 


.100 


.100 


.200 


1.00 


2.00 


Maximum 

Thermal 

Resistance 

(°C/watt) 


60 


40 


60 


30 


30 


20 


20 


Power (max) 
Dissipation 


2.5 


3.75 


2.5 


5.0 


5.0 


7.5 


7.5 



at 25 °C 

(watt) 
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IV. EVALUATION OF SEMICONDUCTOR MICROWAVE SWITCH 

CHARACTERISTICS 



In this chapter, semiconductor device figures of merit, Kawakami’s M and diode 
Q, are tested for the shunted switch configuration, to find the best parameter for use in 
a switch selection and/or design problem. From equation 1.5, we know that the figure 
of merit M, is a function of insertion loss, S 2 1), in the switch on-state, and load 
isolation, S 2 ^( 2 ), in the switch off-state. In order to solve the equation mentioned 
above, the switching characteristics, insertion loss and load isolation, should be 
provided. By using a CAD program (Touchstone) we can minimize the insertion loss in 
the switch on-state, and maximize the load isolation in its off-state. From the value of 
these two switching characteristics, we can compare and select the figure of merit that 
best characterizes the performance of this type of microwave switch. The pow’er 
dissipation and heat removal requirements in the individual semiconductor device can 
be estimated from its terminal voltage. We will introduce two methods to evaluate the 
semiconductor device terminal voltage, and its power dissipation: 

1. Using a Touchstone application program to evaluate the terminal voltage of 
interest. 

2. Using a Fortran program written from the circuit model and method introduced 
in Chapter II. 

A. INSERTION LOSS AND LOAD ISOLATION CALCULATION USING 
COMPUTER-AIDED CIRCUIT ANALYSIS PROGRAM (TOUCHSTONE) 

The Touchstone program is advanced software for RF/ microwave computer- 

aided engineering [Ref. 11]. It can provide an output of the S-parameters of a multi- 

port network and can perform an interactive optimization, and has applications which 

can be used in the measurement of selected terminal voltages. For a given electrical 

circuit, we can create a circuit file and evaluate circuit performance or use the 

optimization mode to get the optimized circuit-component values. 

There are two alternative methods of approach to analyze these two switching 

characteristics. The first method uses semiconductor device parameters, diode junction 

resistance (Rj) and diode junction capacitance (Cp, as the direct input parameters of 

the circuit file. The second method employs semiconductor device figure of merit, diode 

Q, as an input parameter of the circuit file. In this approach the diode junction 
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capacitance and/or resistance are calculated in the circuit file, and then follow the same 
procedure as the first category. 

Two types of microwave switches, single pole single throw, and single pole 
double throw switches, with one or multiple-semiconductor devices are investigated. By 
using the electrical equivalent circuit model introduced in Chapter II, we can create the 
circuit file following one of the above categories for the Touchstone system. The 
assumptions for the circuit model are: 

1. All diodes in the switch are identical. 

2. The transmission line configurations in the switch are symmetrical.' 

3. The transmission line sections are lossless. 

4. The transmission line impedance is between 30 to 90 Ohms. 

5. The transmission line electrical length is between 30 to 120 degrees. 

If the input of the circuit file is the figure of merit, or diode Q, then the junction 
resistance and junction capacitance are calculated from equations 4.2 and 4.3: 



Qdiode 



1000 

2*F 0 RjCj 



1000 

27tFQCjQ c j- oc j e 



(eqn 4.1) 



(eqn 4.2) 



1000 

27tF()CjQ c jj oc j e 



(eqn 4.3) 



where is the figure of merit. 

Rj is the diode junction resistance in ohm. 

Cj is the diode junction capacitance in pF. 

Fq is the center frequency in GHz. 

The relation of Rj, C ; and diode Q at the center frequency of 15 GHz. from equation 
4.2, 4.3 is shown in Figure 4.1. in our circuit file for convenience we fix the value of 
capacitance and calculate the value of junction resistance (equation 4.2) and vary' the 
value of diode Q. 



7 For two-diode switches, TLIN1 is identical to TLIN3 (see Figures 4.3 and 4.7). 
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RESISTANCE VS CAPACITANCE 




Figure 4.1 Resistance vs Capacitance as a Function of Diode Q. 
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FREQUENCY 15 GHZ 



In the procedure of creating a circuit file to optimize the switch characteristics 
for low insertion loss and high isolation, one more assumption is made; the impedance 
values of the diode are kept unchanged at the particular frequency in each state. The 
length and the characteristic impedance of the transmission line sections are varied to 
achieve the optimum switching characteristics in the operating frequency interval. In 
the optimization mode of the circuit file, we can minimize the switch insertion loss and 
maximize load isolation at the same time. This approach may require a long 
computation time to produce a result, and the switch insertion loss still may not 
approach the wanted value. It is found that a more efficient method of approach is to 
optimize only the insertion loss. The computation time is typically shorter, and 
normally produces satisfactory insertion -loss and load isolation. Figure 4.2 is the flow 
chart of steps in creating a circuit file for the Touchstone system. 

1. Insertion Loss in the Switch On-state 

In the switch on-state, the ratio of the available power to the power reaching 
the load is defined as the insertion loss of the switch. In this switching state, the p-i-n 
diode is equivalent to a very high impedance in shunt with the circuit. A small amount 
of energy will be dissipated in the diode, but most of the energy is delivered to the load. 
Therefore the insertion loss of a well designed switch should be as low as possible. In 
order to achieve the low insertion loss, the optimization mode in the Touchstone 
system is used to perform the minimization of insertion loss to obtain the circuit design 
values. In practice the insertion loss has units of dB, and can be calculated from: 



Insertion loss in dB 



101og 10 



[Pjnl 

^load^ 



(eqn 4.4) 



From Chapter II, under the assumption of an ideal switch with matched load at both 
input and output ports, we can use the scattering matrices to represent the 
transmission behavior for both types of switches. Recall from Chapter II, for the on- 

state representation of S-matrices: 



SPST ON STATE: [S] 



0 1 
1 0 



(eqn 4.5) 
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SPDT ARM (2) ON: [S] = 



(eqn 4.6) 



0 1 0 
1 0 0 
0 0 0 



From equations 4.5 and 4.6, the insertion loss of the switch is defined as |S^| in dB 8 
and can be evaluated from the circuit file of Touchstone system. 

2. Load Isolation in the Switch Off-state 

In the switch ofF-state, the load isolation of the switch is also defined by the 
ratio of the available power to the power reaching the load. In this switching state, the 
p-i-n diode is equivalent to a very low impedance in shunt with the circuit. Most of the 
energy is reflected back to the source, some of the energy is dissipated in the diode and 
a small amount of energy can transfer to the load. For an ideal switch in its off-state, 
the load isolation should be infinite, meaning that no energy is transferred to the load. 
In the real switch, the higher load isolation is the better switch performance. The load 
isolation has units of dB and has the same formula as in equation 4.4. In order to 
achieve the best performance of the designed switch, the Touchstone optimization 
mode is used, if desired, to optimize the load isolation as well as insertion loss. The S- 
matrix representation of the switch in the off-state with the same assumptions as for 
the insertion loss is defined as follows: 



SPST OFF STATE: [S] 



1 0 
0 1 



(eqn 4.7) 



SPDT ARM (3) ON: [S] = 



0 0 1 
0 0 0 
1 0 0 



(eqn 4.S) 



From [Refs. 5.6], and equations 4.7 and 4.S |S->j| in dB is denned as the load isolation 
of the SPST switches, and ] S- j | in dB is defined as the load isolation of the SPDT 

switches. 



8 dB value = 20 log jq S 21 
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Figure 4.2 Flow Chart of a Circuit File. 
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3. Example of Insertion Loss and Load Isolation Evaluation 

We divide the evaluation of switching characteristics into two cases: 

1. Single pole single throw switch with two diodes. 

2. Single pole double throw switch with four diodes. 9 

Two categories of creating the circuit file are applied to both cases. 
a. SPST Two-Diode Switch 

The electrical circuit representation of the SPST two-diode switch is showm 
in Figure 4.3(a). The switch on- and off-state equivalent circuits are shown in Figures 
4.3 (b) and (c) respectively. To create a circuit file from the equivalent circuit, node 
numbers are assigned to every component, then the steps in Figure 4.2 are followed. 
The example of circuit files for switching characteristics, insertion loss and load 
isolation, evaluation of a SPST two-diode switch are listed in Appendix A. 

SPSTD.CKT, listed in Appendix Al, is an example of the circuit file 
created for the single pole single throw switch in the first category, in which diode 
resistance and capacitance values are chosen and used as input data for the circuit file. 
Input parameters are diode junction resistance and diode junction capacitance from 
typical p-i-n diode parameters available in Table 1. We select the diode model number 
MA-4P203 with 0.15 pF junction capacitance and 1.5 ohms junction resistance to be 
the input parameters for the circuit file. After optimizing the circuit file from 10 to 20 
GHz, we have the output file of the switch insertion loss and load isolation as showm 
in Table 2. 

From the Touchstone output file listing in Table 2, LINDON and 
LINDOFF are swatch insertion loss and load isolation respectively. Our required 
switch operating frequency range is between 10 and 20 GHz, the wanted value of 
switch insertion loss is around -0.1 to -0.2 dB and the load isolation greater than -45 
dB. The characteristics of Table 2 satisfy our switch design requirement. The plot of 
insertion loss and load isolation is shown in Figure 4.4. 

After optimization, the circuit-component values, characteristic impedances 
and electrical lengths of the transmission line sections, are given in the optimized 
circuit file. The optimized line-section parameters of the example in Appendix Al are: 
transmission line section A has Z a of 39.40 ohms and E a of 106.3 degrees, transmission 
line section B has Z a of 38.25 ohms and of 62.7 degrees. In switch construction, to 
obtain the design value of the insertion loss and load isolation, the transmission line 



9 There are two diodes in each branch of the switch. 
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4.3(a) SPST 2-DIODE EQUIVALENT CIRCUIT 
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4.3(b) SPST ON-STATE EQUIVALENT CIRCUIT 




(C) SPST 0 F 7 -STATE EQUIVALENT CIRCUIT 



Figure 4.3 Single Pole Single Throw Switch Equivalent Circuit. 
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Figure 4.4 SPST Switch Insertion Loss and Isolation vs Frequency. 
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TABLE 2 

TOUCHSTONE OUTPUT SPSTD.CKT 



Frequency - GHz 


dB [S21] 
LINDON 


dB [S21] 
LINDOFF 


10.0000 


-0.150 


-46.276 


12.0000 


-0.053 


-47.025 


14.0000 


-0.0S5 


-48.011 


16.0000 


-0.110 


-49.261 


1S.0000 


-0.141 


-50.590 


20.0000 


-0.195 


-51.750 



sections in the switch should have their characteristic impedances and electrical lengths 
close to the calculated values. 

SPSTQ.CKT, listed in Appendix A2, is an example of a circuit file created 
for the single pole single throw switch in which input parameters are Q^iode anc * diode 
junction capacitance. The selected value of junction capacitance is 0.15 pF and the 
value of diode Q are between 10 and 500 are used as input parameters of this circuit 
file. After optimizing the circuit file from 10 to 20 GHz we have the output file of the 
switch insertion loss and load isolation as shown in Table 3. 

Table 3 is the collection of of SPSTQ.CKT outputs evaluated for different 
diode Q values. This table verifies that the diode Q has the character of a figure of 
merit, serving as an indicator of switching performance. The higher diode Q is, the 
better the switch characteristics: low insertion loss and high isolation. 

Figures 4.5 and 4.6 show tne plot of insertion loss vs frequency and load 
isolation vs frequency of a SPST switch with two diodes with diode Q as a parameter. 
We can use these two figures together with Figure 4.1 for a design problem of a SPST 
switch from a given specification. For example, a specification may require a two p-i-n 
diode SPST switch to have at least -0.1 dB insertion loss and isolation better than 
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TABLE 3 

TOUCHSTONE OUTPUT SPSTQ.CKT 
INSERTION LOSS IN dB 



Frequency 

GHz 


Diode 
Q= 10 


Diode 
Q= 20 


Diode 
Q= 50 


Diode 
Q= 100 


Diode 
Q = 200 


Diode 
Q= 500 


10.0000 


-.367 


-.342 


-.135 


-.064 


-.056 


-.052 


12.0000 


-.194 


-.111 


-.049 


-.031 


-.019 


-.011 


14.0000 


-.342 


-.170 


-.082 


-.062 


-.044 


-.034 


16.0000 


-.499 


-.231 


-.101 


-.052 


-.037 


-.013 


18.0000 


-.517 


-.269 


-.134 


-.072 


-.039 


-.019 


20.0000 


-.724 


-.393 


-.186 


-.096 


-.048 


-.019 


Frequency 

GHz 


Diode 
Q= 10 


LOAD 
Diode 
Q= 20 


ISOLATION IN dB 
Diode Diode 

Q = 50 Q= 100 


Diode 
Q = 200 


Diode 
Q = 500 


10.0000 


-20.21 


-30.53 


-47.38 


-58.78 


- 70.62 


-86.42 


12.0000 


-20.51 


-30.84 


-48.06 


-59.46 


-71.30 


-87.09 


14.0000 


-21.05 


-31.37 


-48.88 


-60.53 


-72.37 


-88.17 


16.0000 


-21.89 


-32.27 


-50.17 


-61.90 


-73.76 


-89.57 


18.0000 


-22.93 


-33.42 


-51.39 


-63.26 


-75.14 


-90.96 


20.0000 


-23.95 


-34.57 


-52.40 


-64.28 


-76.18 


-92.01 



-60.00 dB at the mid-frequency of 15 GHz. In Figure 4.5, with the value of diode Q 
50 or more the switch can have insertion loss less than -0.1 dB. But from Figure 4.6 
the diode Q of 50 is not great enough to obtain load isolation of -60.00 dB. 
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DIODE Q VS INSERTION LOSS 




Figure 4.5 SPST Switch Insertion Loss vs Q< 3 i 0 d e * 
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SPST 2-DIODIC SWITCH 



DIODE Q VS LOAD ISOLATION 
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Figure 4.6 SPST Switch Load Isolation vs Q^q^q- 
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SPST 2-DIODK SWITCH 



Therefore we select a diode for the switch that has the value of Q^j 0( j e of 100. 
From Figure 4.1 with the Q= 100 curve we can choose one of the p-i-n diodes that 
has the value of junction capacitance and resistance close to our estimated value. If 
we select p-i-n diode model number MA-4P504 with diode junction capacitance of 0.2 
pF and junction resistance of 0.6 ohm, at 15 GHz we obtain the diode Q vaiue of 
88.42, switch insertion loss of -.102 dB and isolation of -60 dB from the modified 
SPST.CKT circuit file. 

b. SPDT Four-Diode Switch 

A single pole double throw switch is considered as two SPST.switches with 
input arms tied together, when one arm is in the on-state, the other is in the off-state. 
A SPDT four-diode switch has the electrical circuit representation and the switch on 
off-state equivalent circuit shown in Figure 4.7(a) and (b) respectively. To create a 
circuit file from the equivalent circuit, the same procedure can be followed as in SPST 
switch case. Or we can make some necessary modifications to the SPST switch circuit 
files in Appendix A. The example of circuit files for evaluating SPDT switch insertion 
loss and load isolation are listed in Appendix B. 

SPDTD.CKT, listed in Appendix Bl, is an example of the circuit file used 
for the single pole double throw' switch with four diodes. With the same p-i-n diode 
model number as in the circuit file SPSTD.CKT, after optimizing this circuit file from 
10 to 20 GHz, we have the output file of the insertion loss and load isolation shown in 
Table 4. The optimized transmission line section characteristic impedances and 
electrical lengths are: 46.56 ohms and 93.58 degrees for transmission line section A, 
89.06 ohms and 32.82 degrees for transmission line section B. 

The output file listing in Table 4 shows that for the required operating 
frequency between 10 to 20 GHz, the designed SPDT four-diode switch using p-i-n 
diode model number MA-4P203 with 0.15 pF junction capacitance and 1.5 ohms 
junction resistance, has insertion loss around -1.0 dB and load isolation greater than 
-50 dB We realize that the insertion loss of the SPDT switch case is not so low as for 
the SPST switch case. How'ever the insertion loss and load isolation satisfy the general 
requirement of a switch design problem. The piot of insertion loss and load isolation 
vs operating frequency is shown in Figure 4.8. 

SPDTQ.CKT, listed in Appendix B2, is an example of the circuit file 
created for the single pole double throw' sw’itch for which input parameters are 
and diode junction capacitance. In this example, the value of 0.15 pF junction 
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Figure 4.7 Single Pole Double Throw Switch Equivalent Circuit. 
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Figure 4.S SPDT Switch Insertion Loss and Isolation vs Frequency. 
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TABLE 4 

TOUCHSTONE OUTPUT SPDTD.CKT 


Frequency - GHz 


dB [S21] 


dB [S21] 




LINDON 


LINDOFF 


10.0000 


-0.756 


-56.199 


12.0000 


- 0.327 


-57.856 


14.0000 


-0.226 


-59.898 


16.0000 


-0.385 


-61.603 


1S.0000 


-0.417 


-61.471 


20.0000 


-0.665 


-59.573 



capacitance, and values of diode Q from 10 to 500 are selected and used as input 
parameters of the circuit file. The optimized output file SPDTQ.CKT is shown in 
Table 5. 

For general use of diode Q as the figure of merit to characterize the 
switching performance of this type of switch, the choice of junction capacitance value 
selected as input parameter to the Touchstone file does not strongly affect the Q- 
dependence of switch performance shown in Table 5. The result of output file 
SPDTQ.CKT listed in Table 5 verifies that the diode Q has the character of a figure of 
merit, serving as an indicator of the switch performance. In a switch design problem, 
diode Q can be used as a parameter that characterizes the switching transmission 
characteristics. The plot of insertion loss and load isolation vs operating frequency as 
a function of diode Q from the results in Table 5 are shown in Figures 4.9 and 4.10. 
These two piots can be used in a design problem of a SPDT switch from a given 
specification. The procedures for the design of a SPDT switch are the same as 
mentioned in the SPST switch case. 



51 



TABLE 5 

TOUCHSTONE OUTPUT SPDTQ.CKT 



INSERTION LOSS IN dB 



Frequency 

GHz 


Diode 
Q= 10 


Diode 
Q= 20 


Diode 
Q = 50 


Diode 
Q= 100 


Diode 
Q = 200 


Diode 
Q= 500 


10.0000 


-1.96 


-1.38 


-1.04 


-.926 


-.824 


-.796 


12.0000 


-1.36 


-.725 


-.345 


-.217 


-.148 


-.117 


14.0000 


-1.13 


-.577 


-.265 


-.179 


-.129 


-.095 


16.0000 


-1.32 


-.756 


-.394 


-.294 


-.232 


-.190 


1S.0000 


-1.69 


-.935 


-.395 


-.214 


-.118 


-.059 


20.0000 


- 2.37 


-1.49 


-.835 


-.564 


-.515 


-.425 


Frequency 

GHz 


Diode 
Q= 10 


LOAD ISOLATION 
Diode Diode 

Q= 20 Q= 50 


IN dB 
Diode 
Q= 100 


Diode 
Q = 200 


Diode 
Q = 500 


10.0000 


-33.55 


-44.74 


-59.50 


-71.18 


-83.59 


-99.46 


12.0000 


-34.35 


-45.92 


-60.83 


-72.58 


-85.09 


-100.95 


14.0000 


-35.82 


-47.90 


-63.08 


-74.99 


-87.56 


-103.43 


16.0000 


-37.56 


-49.71 


-64.75 


-76.66 


-89.13 


-105.00 


18.0000 


-38.27 


-49.68 


-64.27 


-75.97 


-88.18 


-105.04 


20.0000 


-37.54 


-48.24 


-62.75 


-74.37 


-86.21 


-102.20 



B. THE SEMICONDUCTOR DEVICE FIGURE OF MERIT WHICH BEST 
CHARACTERIZES THE SWITCHING PERFORMANCE 

Two well known figures of merit were used as parameters to predict the switching 

performance of semiconductor microwave switches. In this thesis research, the 
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Figure 4.9 SPDT Switch Insertion Loss vs Q^iode' 
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SPDT 4-DIODE SWITCH 
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Figure 4.10 SPDT Switch Load Isolation vs Q<^ oc j e - 
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SPDT 4— DIODE SWITCH 



semiconductor device considered is the p-i-n diode, used in a shunt-configuration 
microwave switch. The conventional diode Q was shown to coincide with the 
simplified equation of Kurokawa and Schlosser for quality factor Q, and therefore we 
use the conventional Q as our semiconductor device figure of merit. In Chapter II 
equations 2.10 and 2.11 showed the relationship between both figures of merit: the 
diode Q and Kawakami’s M. To verify these two equations, we investigated the 
relationship of Kawakami’s M and the switch transmission characteristics defined in 
equation 1.5. Optimal values of switch insertion loss and load isolation were evaluated 
from the Touchstone circuit file which has the diode Q as an input parameter. For the 
SPDT switch case, values of insertion loss and load isolation evaluated with different 
diode Q values at 15 GHz. are substituted in equation 1.5, and the results are listed in 
Table 6. 



TABLE 6 

COMPARISON BETWEEN FIGURES OF MERIT M, AND DIODE Q 



Diode Q 


S 21 (1) 


S 2 i (2) 


M (eqn 1.5) 


AM 

— x 100 (%) 
A Q 


10.0000 


0.9421413 


0.1590319 


0.9770399 




20.0000 


0.9677804 


0.0865225 


0.9937131 


0.1667 


50.0000 


0.9830945 


0.0404859 


0.9985498 


0.0161 


100.000 


0.9872322 


0.0222759 


0.9994510 


1.8e-3 


200.000 


0.9897033 


0.0118343 


0.9997507 


2.9e-4 


500.000 


0.9915362 


0.005379S 


0.9999076 


5.2e-5 



Table 6 shows the comparison of the two figures of merit, with M calculated 
from equation 1.5. Note that in this table S 21 O), S2j(2) correspond to switch 
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insertion loss and load isolation respectively, where the value of S 21 is in the 
form: S 2 j = lO^k/20). The values of M are less than 1 as mentioned earlier. For 

this type of switch at Qj- oc j e of 10 the value of M is 0.9762, and at of 500 the 

value of M is 0.9999. For typical diode Q values, M departs very little from unity, and 
the percentage of deviation is reduced tremendously in the high Q region. Thus the 
parameter is more convenient to use due to its wider range compared with the 

parameter M. When tested for different types of switch, 10 the results are nearly the 
same as shown in Table 6 for all the values of M (almost unity and small deviation in 
each interval). An additional reason in support of our decision to select the 
conventional diode Q as the best figure of merit to characterize switching performance 
is that it is well knowm and widely used. Table 7 shows the comparison of the result 
from Table 6 and values of M calculated from equations 2.10 and 2.11, this table 
verifies that equation 2.10 can be used to convert the diode Q value to the figure of 
merit M, and vice versa. 



TABLE 7 

COMPARISON OF FIGURES OF MERIT M, FROM DIFFERENT 

APPROACHES 



Diode Q 


M (eqn 1.5) 


M (eqn 2.10) 


M (eqn 2.11) 


10.0000 


0.9770399 


0.9805806 


0.9800000 


20.0000 


0.9937131 


0.9950372 


0.9950000 


50.0000 


0.9985498 


0.9992001 


0.9992000 


100.0000 


0.9994108 


0.9998000 


0.9998000 


200.0000 


0.9997507 


0.9999499 


0.9999500 


500.0000 


0.9999076 


0.9999919 


0.9999920 



10 SPDT switch with tuning stub, microstrip model. 
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C. NODAL VOLTAGE CALCULATION USING TOUCHSTONE PROGRAM 



The power dissipation in a particular diode can be calculated by the voltage 
across its resistance or the current flow in the resistor of the p-i-n diode equivalent 
circuit. In the Touchstone program there is a technique for probing the nodes within a 
circuit file to determine the voltages at those nodes. An ideal voltage-controlled 
voltage source (VCVS) is used as a non-loading voltmeter probe [Ref 12]. The input 
terminals of the VCVS are connected to the nodes between which the voltage is to be 
determined. The output terminals are connected to an output port so that the voltages 
may be read. The voltages read out from this method are normalized to unity value of 
the open circuit generator voltage [Ref. 12]. Appendix C shows examples of the 
Touchstone circuit files written for the evaluation of nodal voltages of SPST and SPDT 
switches. For SPDT switch case, the equivalent circuit model using in nodal voltage 
calculation is shown in Figure 4.11. 

The circuit file SPDTV.CKT, listed in Appendix C2, is an example of a 
Touchstone circuit file using VCVS modules acting as voltage probes reading the nodal 
voltages from the attached terminals. The input parameters for this circuit file are the 
diode parameters and the optimal transmission line section values of the SPDT switch 
evaluated in Appendix B2. The major disadvantage of the Touchstone circuit file is 
that we can define up to four ports in our circuit file, meaning that we can measure no 
more than three nodes at a time. To measure a SPDT switch with four diodes, there 
are six nodes of interest. The first three nodes are measured in the SPDTV.CKT circuit 
file, while the last three nodes are measured by making some modifications to this 
circuit file (by moving the VCVS module to the new positions). The output files of 
SPDTV.CKT are listed in Table 8. 

Knowing terminal voltage across the resistive part of the p-i-n diode equivalent 
circuit both on and off-state, the next step is to calculate the power dissipation in the 
individual p-i-n diode from the equation: 







(eqn 4.9) 



where P n is the dissipated power at node n, 

V n is the terminal voltage at node n, and 

R n is the resistive component of the n^ p-i-n diode. 
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Figure 4.11 Nodal Voltage Measurement Equivalent Circuit. 



TABLE 3 

NODAL VOLTAGE FROM TOUCHSTONE SPDTV.CKT 

SWITCH ON-STATE SWITCH OFF-STATE 



Frequency 


"^load 


V D1 


V D2 


^load 


V D1 


V D2 


10.0000 


.458 


.006 


.006 


7.7e-4 


.017 


7.3e-4 


12.0000 


.482 


.008 


.008 


6.4e-4 


.016 


6.0e-4 


14.0000 


.4S5 


.010 


.009 


5. le-4 


.015 


4.7e-4 


16.0000 


.478 


.012 


.010 


4.2e-4 


.014 


3.9e-4 


18.0000 


.477 


.015 


.011 


4.2e-4 


.015 


4.0e — \ 


20.0000 


.463 


.018 


.012 


5.3e-4 


.021 


5.0e-4 



5S 



In practice the dissipated power in a particular diode is normalized to the power 
delivered to the load in the switch on-state. Therefore each power dissipation calculated 
from equation 4.9 is normalized to the switch on-state output power. For the 
normalized power value, equation 4.9 becomes: 

2 

p _ Vfl R joad . (eqn 4.10) 

n ivS 

where Vj oa( j is the load terminal voltage of the switch on-state. 

The use of Touchstone circuit files, SPSTV.CKT and SPDTV.CKT, is very 
convenient for evaluating the diode terminal voltage, but we cannot use Touchstone to 
solve for power dissipation and the diode junction temperature problem directly. To do 
so one may use a hand calculator or write a special computer program calculating the 
power dissipation in each diode from equation 4.10, and use these dissipated powers to 
calculate the heat production from equations in Chapter III. However we see that the 
Touchstone output file has only three significant digits after the decimal point, 
meaning that numerical round off takes place in the third digit. These outputs 
generally satisfy engineering requirements. However for a very small nodal voltage 
across the diode, these rounded off outputs may cause significant error when we 
substitute them in a quadratic equation such as equation 4.10. It is a serious 
disadvantage to use these outputs for high output power microwave switches. For 
example an output voltage of 0.010 from Touchstone to represent a value of 0.0103 
volts can cause an error of 0.88 watts of diode power in a one kilowatt switch. This 
number may look small, however it can cause an error of 30-40 °C in diode junction 
temperature rise. This error can affect the decision of selection of a proper p-i-n diode 
for the switch design. 

D. EVALUATION OF POWER DISSIPATION AND HEATING USING 

FORTRAN PROGRAM 

To avoid some disadvantages of the use of the Touchstone nodal voltage 
evaluation for power computation, we may search for another CAD program or create 

a specific computer program to solve this problem. Many microw'ave CAD programs 
have been written in FORTRAN IV for both mainframe and minicomputers, but up to 
now there is no program written to solve for the power dissipation in a microwave 
circuit [Ref. 7]. The Fortran language requires specific formats of input and output 
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and has to be compiled before running the program. Nonetheless, the advantages to 
the microwave engineer are its complex number notation, and double precision 
capability. These two advantages lead to accurate outputs, for example in the 
multiplication of more than ten complex matrices in cascade. To write a Fortran 
program to evaluate the terminal voltages of the diode microwave switch equivalent 
circuit, electrical two-port components of the switch are converted _ into simple 
appropriate two-port matrix forms. In chapter 2 we introduced an ABCD matrix 
representation of the two-port networks for a transmission line section, a shunt 
admittance and series impedance in the lumped-circuits. For our switch circuit model, 
the circuit components are characterized as boxes in cascade connecting with the 
applied voltage at one end and the load at the other end. 

To evaluate the power dissipation and junction temperature in the diodes of a 
SPDT switch which has an equivalent circuit shown in Figure 4.12, with the optimal 
transmission line sections by use of Touchstone, the procedure follows the flow 
diagram shown in Figure 4.13. 

1. Nodal Voltage and Current Calculation 

We assume that the switch is supplied by a one volt generator with 50 ohms 
series impedance. First we divide the equivalent circuit into two parts, one the on- 
state branch the other the off-state branch. For both on- and off-states, each 
transmission line section and the diode equivalent circuit two-port network is 
transformed into ABCD matrix form (Figure 4.12(a)) by subroutine MATRIX. This 
subroutine using equations 2. IS and 2.19 for the transformation of a shunt diode and a 
lossless transmission line section into ABCD two-port matrices respectively. To obtain 
the voltages and currents in the circuit, the successive input impedances of the 
individual two-port elements must be determined. All two-port networks in the branch 
are cascaded by subroutine MATMUL and produce a new equivalent circuit as in 
Figure 4.12(a). Subroutine MATMUL provides a complex matrix product of a pair of 
matrices. From Figure 4.12(b), input impedance in each branch (Zj non , Z- 0 fy) is 
calculated from equation 2.21. The impedance of the switch (Z- n ) is determined by the 
paralleling of input impedances of both branches lZ; non , Z; no £f). The switch input 
current (Ij n ) and input voltage (v- n ) from the generator are calculated from equations 
2.22 and 2.23. The next step is to find the currents and voltages in the on- and off- 
state branches. The input voltage on each branch is equal to the switch input voltage 
(paralleled ports). Current flow in the on- and off-state branches is determined by 
current divider method as follows: 
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Figure 4.12 ABCD Matrix Representation of SPDT switch Circuit Model. 
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OFF-STATE BRANCH 
4.12(b) 






Figure 4.13 Flow Diagram for Calculation of Power Dissipation 
and Diode Junction Temperature. 
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on 




(eqn 4.1 1) 





(eqn 4.12) 




where I Qn is the current flow in switch on- state branch, 

I is the current flow in the switch ofF-state branch, 

Z^on * s i n P ut impedance of the switch on-state branch and, 

Z^ff is the input impedance of the switch off-state branch. 

The final step of nodal voltage and current calculation returns to Figure 
4.12(a). In each branch, starting from the leftmost node, currents and voltages in 
adjacent nodes are calculated by equations 2.24 and 2.25. In this step we obtain the 
currents flowing in and out of each particular node, as well as the input and output 
voltages, leading us to the power calculation as in the next section. 

2. Power Dissipation in the Individual P-I-N Diode 

For the calculation of diode power dissipation in each branch of the switch, 
we create the subroutine named POWER. In this subroutine, the input currents and 
input voltage of each branch are the initial input parameters. These currents together 
with input and output voltages of the diode’s node are calculated from equations 2.24 
and 2.25. Equation 2.26 is used to calculate the diode power dissipation with its 
corresponding nodal voltage and current. For the pow r er at load we can use the 
rightmost nodal current and voltage (I n , V n ) as inputs to equation 2.27 or using 
equation 2.2S with the branch input power minus the total of all power dissipation in 
that branch (include power dissipation in the transmission line sections). The power 
dissipation in each diode and all the calculated powers are normalized to the power 
delivered to the load in the switch on-state. 

3. Junction Temperature in the Individual P-I-N Diode 

The junction temperature of a p-i-n diode is determined by the ambient 
temperature of the switch and the power dissipation in the diode itself. The ambient 
temperature is assumed to be constant at about the room temperature. The power 
dissipation in the diode is a major factor that causes damage to the semiconductor 
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devices in the switch. From the calculation of power dissipation mentioned in previous 
section, we can separate the calculation of diode junction temperature into two cases, 
CW power dissipation and pulsed power dissipation. In the Fortran program, the 
subroutine THERM 1 evaluates the junction temperature of the CW-power case, for 
the diode mounted on a heat sink. The junction temperature is considered by summing 
the ambient temperature .and the increment of temperature of the diode caused by the 
input CW-power (equation 3.7). The total thermal resistance, 0j c , of a particular 
diode is available from the p-i-n diode specification sheet. For the pulsed input-power 
case, the subroutine THERM2 is used for evaluating the diode junction temperature 
rise. In this case we assume that the input pulse width is small compared with the 
thermal time constant, meaning that all the thermal energy is absorbed in the diode, 
with no energy spreading out to the surroundings. The junction temperature is 
obtained by summing the ambient temperature and the temperature rise of the junction 
(equation 3.12), where the temperature rise of the junction is determine by the thermal 
capacity of the diode active region (HC), the pulse width in microsecond and the power 
dissipation in the diode (equation 3.11). For the switch design problem, the 
recommended maximum junction temperature for safe operation is 200 °C 
[Ref. 10: p.290]. 

4. An Example and Results of Power Dissipation and Diode Junction Temperature 

Evaluation 

Appendix D is the listing of two Fortran programs using for the evaluation of 
power dissipation and junction temperature in each diode of a single pole single throw 
switch, and a single pole double throw switch with one or more diodes. These two 
programs, PSPST and PSPDT, are written for general use. The number of diodes in 
each branch is an input variable that determines the switch construction (number of 
transmission line sections). Diode input parameters can be either diode Q or junction 
resistance and junction capacitance. All input parameters can be entered interactively 
or in an unformatted input file, and the results of the evaluation are listed in an output 
file. As an example of power dissipation and junction temperature in a SPDT.switch 
with four MA--1P203 p-i-n diodes, we use PSPDT program in Appendix D2 for the 
evaluation. From the specification sheet, diode MA-4P203 has 0.15 pF. junction 
capacitance, 1.5 ohms junction resistance and 30° C/watt total thermal resistance. The 
value of transmission characteristic impedance and electrical length for each 
transmission line section are the optimal component values obtained from 
SPDTD.CKT Touchstone circuit file. For a 100 watt CW power SPDT switch 
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operating at 10 to 20 GHz frequency range, the results of using PSPDT program are 
listed in Table 9. 



TABLE 9 



DIODE POWER DISSIPATION AND JUNCTION TEMPERATURE OF 
THE 100- WATT CW POWER SWITCH 

ON-STATE OFF-STATE 



FRE 


POWER(watts) 


TEMP(°C) 


POWER(watts) 


TEMP(°C) 


(GHz) 


D 1 


d 2 


D 1 


d 2 


D 1 


d 2 


D 1 


d 2 


10.00 


0.814 


0.578 


39.4 


32.3 


2.523 


0.002 


90.7 


15.1 


12.00 


1.124 


0.832 


48.7 


40.0 


2.653 


0.002 


94.6 


15.1 


14.00 


1.486 


1.133 


59.6 


48.9 


2.770 


0.002 


98.1 


15.1 

+ 


16.00 


1.908 


1.479 


72.2 


59.4 


2.865 


0.002 


100.9 


15.1 


18.00 


2.407 


1.872 


87.2 


71.2 


2.932 


0.002 


102.9 


15.1 


20.00 


3.003 


2.311 


104.1 


84.3 


2.971 


0.002 


104.1 


15.1 



From the result in Table 9, we realize that most of the switch power is 
dissipated in the first diode of each branch, particularly in the switch-off branch. In the 
switch on-state the power is shared among the two diodes, with the first one dissipating 
more power. For the whole range of frequencies the first diode of the off-state branch 
dominates the power dissipation and undergoes the greatest junction temperature rise. 
The switch can be damaged if the junction temperature exceed its maximum allowable 
value (200 °C ). In most diode specification sheets, the maximum power dissipation at 
the given ambient temperature is given based on the CW-input power to the diode. 
The CW-input power has more effect in the temperature rise in the diode junction 
because of the heat accumulated continuously in the junction. For this reason a diode 
switch with short-duration input power pulses can handle higher peak powers than for 
CW-input power. For example the MA-4P203 diode has the maximum power handling 
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capability of 5 watts of CW-input power at 25 °C ambient temperature. In Table 9 the 
designed switch has a maximum power dissipation in the diode of about three watts 
and 105°C junction temperature (20°C ambient temperature). The result of Table 9 
shows that our switch design is capable to handle a 100-watt input CW power. Figure 
4.14 shows the plot of normalized powers to one watt output power of an individual 
diode in both branches. 
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Figure 4.14 Power Dissipation in an Individual Diode. 
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V. CONCLUSIONS AND REMARKS 



In microwave switching design, the CAD program (Touchstone) and Fortran 
programs, are used for evaluating different types of microwave diode switches in order 
to obtain their estimated switching characteristics. The selected diode Q as a 
switching-semiconductor figure of merit can be used as the measure of device switching 
effectiveness, as shown in chapter 4. For radar applications, a single pole double throw 
switch with high isolation over a broad operating frequency range is an essential 
requirement for all radar systems. An addition of more diodes to the switch can 
improve the switch isolation, while increasing the insertion loss. Figures 5.1 and 5.2 
are plots of the switch insertion loss and isolation for switches with one or more 
diodes. The insertion loss and load isolation in these two figures are produced from the 
Touchstone circuit file evaluated for the different diode Q values. These two figures 
are very useful to the designer selecting the number of diodes to be used for specified 
switch insertion loss and isolation, based on the diode Q value. An appropriate Q value 
is selected to characterize the behavior of this switch. The Fortran program is used to 
estimate the power handling capability of the designed switch. The resistor in the diode 
equivalent circuit plays an important role in its power dissipation. The value of diode 
Q, from equation 4.1, is a function of both junction capacitance and resistance. One of 
these two values can be set as an arbitrary parameter, meaning that we can select more 
than one diode having the same diode Q value. The problem is to select that p-i-n 
diode which can handle the power dissipation of the switch. One suggestion is to select 
lower diode resistance to avoid too much power dissipation. However this leads to a 
trade-off in increased power dissipation in the on-state branch. An example of this is in 
the previous chapter; the p-i-n diode model number MA-4P203 has the junction 
resistance of 1.5 ohms and diode Q of 47.15 at 15 GHz, and power dissipation in the 
on-state branch of the first diode is very close to that of the first diode in the off-state 
branch (1.6S and 2. 32 watts). In some switch applications, it is not desirable to have 
heat accumulated in the diode junction at all times (both switching state), which may 
cause shorter diode lifetime. If we select a p-i-n diode which has the same diode Q 
value but greater junction resistance, for example 2.0 ohms junction resistance, we 
obtain 3.478 watts power dissipated in the off-state branch and 1.313 watts in the on- 
state branch. The best way to select the diode for the switch is by using the Fortran 
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program to estimate the power dissipation in the first diode of each branch, by fixing 
the diode Q value and varying diode junction resistance and select the satisfactory 
outcome. 

Due to the approximate circuit model used for the diode microwave switch 
equivalent circuit, the power calculation cannot produce a reliable substitute for 
measurements in the laboratory. Nonetheless it can reduce the design cost by providing 
initial guide lines to the designer. On the other hand these assumptions reduce the 
complexities of the circuit analysis and are mainly employed in the test of different 
figures of merit. Improvements for more accurate result can be achieved as follows: 

1. Using more accurate p-i-n diode equivalent circuit model. 

2. Introducing loss in the transmission line sections. 

From the knowledge that the first diode in the off-state branch dominates the 
power dissipation, in a switch circuit design for more efficient switch performance, the 
diodes used in the circuit can be differed (not necessarily identical). The one closest to 
the generator should have higher diode Q value (lower junction resistance) in order to 
dissipate less power, allowing the remaining diodes to share the power dissipation. 
This method can reduce the maximum power dissipation that occurs in the first diode 
of the off-state branch, and improve the power handling capability of the switch. We 
still can use the Touchstone circuit file to evaluate the switch performance by making 
some necessary modifications to it. The switch optimal insertion loss and isolation are 
processed by the optimization mode to obtain the optimal switch circuit component 
values to be used in the power dissipation evaluation with Fortran program. 
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DIODE Q VS INSERTION LOSS 
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Figure 5.1 Switch Insertion Loss vs Diode Q for Multi-diodes Switch. 
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SPDT DIODE SWITCH 



DIODE Q VS LOAD ISOLATION 
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Figure 5.2 Switch Isolation vs Diode Q for Multi-diodes Switch. 
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APPENDIX A 

TOUCHSTONE CIRCUIT FILE FOR SPST SWITCH 



This appendix provides the listing of two circuit files using in the Touchstone 
CAD program for the evaluation switch insertion loss and load isolation for single 
pole single throw two-diode switches. 

Al. SPSTD.CKT evaluate with diode parameters (C;, Rj). 

A2. SPSTQ.CKT evaluate with diode figure of merit (diode Q). 



Al. SPSTD.CKT 



! ! SINGLE POLE SINGLE THROW DIODE SWITCH 
! ! INSERTION LOSS AND LOAD ISOLATION CALCULATION 
! ! INPUT PARAMETERS : 

!! 1. JUNCTION CAPACITANCE, CJ 
! ! 2. JUNCTION RESISTANCE, RJ 
! ! 3. CENTER FREQUENCY, FO 
! .' OUTPUT : 

! ! OPTIMIZED TRANSMISSION LINE IMPEDANCE 
! ! AND ELECTRICAL LENGTH 
! ! MINIMIZED INSERTION LOSS 
! ! LOAD ISOLATION 
VAR 

! INPUT PARAMETERS: 



FO = 


15 


! FREQUENCY IN GHz 


CJ = 


.15 


! JUNCTION CAPACITANCE IN pF 


RJ = 


1.5 


! JUNCTION RESISTANCE IN Ohm 



! OPTIMIZED TRANSMISSION LINE IMPEDANCE (Ohm) 

ZA#30 39.40 90 
ZB#30 38.25 90 

! OPTIMIZED TRANSMISSION LINE ELECTRICAL LENGTH (Degree) 
EA#30 106.3 120 
EB#30 62.70 120 
CKT 

TLIN 1 2 Z A ZA E A EA F A F0 

DEF2P 1 2 SEG1 

TLIN 1 2 Z A ZB E A EB F A F0 

DEF2P 1 2 SEG2 

SRC 1 2 R A RJ C A CJ 

DEF2P 1 2 DON 

RES 1 2 R A RJ 

DEF2P 1 2 DOFF 

! ON STATE CALCULATION 

SEG1 1 2 

DON 2 0 

SEG2 2 3 

DON 3^0 

5EG1 3 4 

DEF2P 1 4 LINDON 

! OFF STATE CALCULATION 

SEG1 1 2 

DOFF 2 0 

SEG2 2 3 

DOFF 3 0 

SEG1 3 4 

DEF2P 1 4 LINDOFF 
OUT 

LINDON DB (S21 ) GR1 ! INSERTION 'LOSS IN dB 
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LINDOFF DB (S21 ) GR1A ! LOAD ISOLATION IN dB 

FREQUENCY 

SWEEP 8 22 2 ! FREQUENCY SWEEP FROM 8 TO 22 GHz 

GRID 

GR1 -2 0 .25 

GR1A -20 -80 

OPT 

RANGE 10 20 

LINDCN DB(S21 )=0 ! MINIMIZE INSERTION LOSS 



A2. SPSTQ.CKT 



! SINGLE POLE SINGLE THROW DIODE SWITCH 
! INSERTION LOSS AND LOAD ISOLATION CALCULATION 
! INPUT PARAMETERS : 

! 1. JUNCTION CAPACITANCE, CJ 
! 2. FIGURE OF MERIT (DIODE Q),Q 
! 3. CENTER FREQUENCY, FO 
! OUTPUT : 

! OPTIMIZED TRANSMISSION LINE IMPEDANCE 
! AND ELECTRICAL LENGTH 
! MINIMIZED INSERTION LOSS 
! LOAD ISOLATION 
VAR 

! INPUT PARAMETERS: 

FO = 15 ! FREQUENCY IN GHz 

CJ = .15 ! JUNCTION CAPACITANCE IN pF 

Q = 100 ! FIGURE OF MERIT (DIODE Q) 

! OPTIMIZED TRANSMISSION LINE IMPEDANCE (Ohm) 

ZA#30 36.66 90 
ZB#30 37.30 90 

! OPTIMIZED TRANSMISSION LINE ELECTRICAL LENGTH (Degree ) 
EA#30 108.3 120 
EB#30 58.7 120 

EQN 

RJ = 1000/ (2*PI*F0*CJ*Q) 

CKT 

TLIN 1 2 Z A ZA E A EA F A F0 
DEF2P 1 2 SEG1 
TLIN 1 2 Z A ZB E A EB F A F0 
DEF2P 1 2 SEG2 
SRC 1 2 R A RJ C A CJ 
DEF2P 1 2 DON 
RES 1 2 R A RJ 
DEF2P 1 2 DOFF 
! ON STATE CALCULATION 
SEG1 1 2 
DON 2 0 
SEG2 2 3 
DON 3 0 
SEG1 3 4 

DEF2P 1 4 LINDON 
! OFF STATE CALCULATION 
SEG1 1 2 
DOFF 2 0 
SEG2 2 3 



! JUNCTION RESISTANCE CALCULATION 



DOFF 3 0 
SEG1 3 4 

DEF2P 1 4 LINDOFF 
OUT 

LINDON DB(S21 ) GR1 
LINDOFF DB(S21) GR1A 
FREQUENCY 
SWEEP 8 22 2 
GRID 

GR1 -2 0 .25 



! INSERTION LOSS IN dB 
! LOAD ISOLATION IN dB 

! FREQUENCY SWEEP FROM 8 TO 22 GHz 



73 



GR1A -20 -80 

OPT 

RANGE 10 20 

LINDON DB (S21 ) = 0 ! MINIMIZE INSERTION LOSS 
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APPENDIX B 

TOUCHSTONE CIRCUIT FILE FOR SPDT SWITCH 



This appendix provides the listing of two circuit flies using in the Touchstone 
CAD program for the evaluation switch insertion loss and load isolation for single 
pole double throw four-diode switches. 



Bl. SPDTD.CKT evaluate with the diode parameters (Cj, Rp. 
B2. SPDTQ.CKT evaluate with diode figure of merit (diode Q). 



Bl. SPDTD.CKT 



! ! SINGLE POLE DOUBLE THROW DIODE SWITCH 
l ! INSERTION LOSS AND LOAD ISOLATION CALCULATION 
! ! INPUT PARAMETERS : 

111. JUNCTION CAPACITANCE, CJ 

112 . JUNCTION RESISTANCE, RJ 
!! 3. CENTER FREQUENCY, FO 

! ! 4 . TERMINAL LOAD 
! ! OUTPUT : 

! ! OPTIMIZED TRANSMISSION LINE IMPEDANCE 
.' ! AND ELECTRICAL LENGTH 
! ! MINIMIZED INSERTION LOSS 
! ! LOAD ISOLATION 
VAR 

1 INPUT PARAMETERS: 

FO = 15 ! FREQUENCY IN GHz 

CJ = .15 ! JUNCTION CAPACITANCE IN pF 

RJ = 1.5 ! JUNCTION RESISTANCE IN Ohm 

RR = 50 ! TERMINAL LOAD IN Ohm 

l OPTIMIZED TRANSMISSION LINE IMPEDANCE (Ohm) 

ZA#30 46.56 90 
ZB#30 39.06 90 

! OPTIMIZED TRANSMISSION LINE ELECTRICAL LENGTH(Deqree ) 
EA#30 93.58 120 
EB#30 32.82 120 
CKT 

TLIN 1 2 Z A ZA E A EA F A F0 

DEF2P 1 2 SEG1 

TLIN 1 2 Z A ZB E A EB F A F0 

DEF2P 1 2 SEG2 

SRC 1 2 R A RJ C A CJ 

DEF2P 1 2 DON 

RES 1 2 R A RJ 

DEF2P 1 2 DOFF 

! ON STATE BRANCH 

SEG1 1 2 

DON 2 0 

SEG2 2 3 

DON 3 0 

SEG1 3 4 

! OFF STATE BRANCH 

SEG1 1 5 

DOFF 5 0 

SEG2 5 6 

DOFF 6 0 

SEG1 6 7 

! DEFINE 3 PORTS SWITCH 
DEF3P 147 SWITCH 
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SWITCH 123 
RES 3 0 R A RR 

DEF2P 1 2 SWON ! 

SWITCH 123 
RES 2 0 R A RR 

DEF2P 1 3 SWOFF ! 

OUT 

SWON DB(S21 ) GR1 ! 

SWOFF DB (521 ) C-R1A ! 

FREQUENCY 

SWEEP 3 22 2 ! 

GRID 

GR1 -2 0 .25 

GR1A -20 -80 

OPT 

RANGE 10 20 

LINDON DB(S21)=0 ! 



INSERTION LOSS SWITCH ON STATE 

LOAD ISOLATION SWITCH OFF STATE 

INSERTION LOSS IN dB 
LOAD ISOLATION IN dB 

FREQUENCY SWEEP FROM 8 TO 22 GHz 



MINIMIZE INSERTION LOSS 



B2. SPDTQ.CKT 



! ! SINGLE POLE DOUBLE THROW DIODE SWITCH 
! ! INSERTION LOSS AND LOAD ISOLATION CALCULATION 
! ! INPUT PARAMETERS : 

! ! 1. JUNCTION CAPACITANCE, CJ 
! ! 2. FIGURE OF MERIT (DIODE Q),Q 
l! 3. CENTER FREOUENCY, FO 
! ! 4 . TERMINAL L5AD , RR 
! ! OUTPUT : 

! ! OPTIMIZED TRANSMISSION LINE IMPEDANCE 
! ! AND ELECTRICAL LENGTH 
! ! MINIMIZED INSERTION LOSS 
! ! LOAD ISOLATION 
VAR 

! INPUT PARAMETERS: 

FO = 15 ! FREQUENCY IN GHz 

CJ = .20 ! JUNCTION CAPACITANCE IN pF 

Q = 100 ! FIGURE OF MERIT (DIODE Q) 

RR = 50 ! TERMINAL LOAD IN Ohm 

! OPTIMIZED TRANSMISSION LINE IMPEDANCE (Ohm) 

ZA#30 36.66 90 
ZB#30 37.30 90 

! OPTIMIZED TRANSMISSION LINE ELECTRICAL LENGTH (Degree) 

EA#30 108.3 120 
EB#30 53.7 120 

EQN 

RJ = 1000/ (2*PI*F0*CJ*Q) ! JUNCTION RESISTANCE CALCULATION 

CKT 

TLIN 1 2 Z A ZA E A EA F A F0 

DEF2P 1 2 SEG1 

TLIN 1 2 Z A ZB E A EB F A F0 

DEF2P 1 2 SEG2 

SRC 1 2 R A RJ C A CJ 

DEF2? 1 2 DOM 

->T7C 1 9 JAOJ 

DEF2? 1 2 DOFF 

! ON STATE BRANCH 

SEG1 1 2 

DON 2 0 

SEG2 2 3 

DON 3 0 

SEG1 3 4 

! OFF STATE BRANCH 
SEG1 1 5 
DOFF 5 0 
SEG2 5 6 
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DOFF 6 0 
SEG1 6 7 

! DEFINE 3 PORTS SWITCH 
DEF3P 147 SWITCH 
SWITCH 123 
RES 2 0 R A RR 

DEF2P 1 2 SWON ! INSERTION LOSS SWITCH ON STATE 

SWITCH 123 
RES 3 0 RARR 

DEF2P 1 3 SWOFF ! LOAD ISOLATION SWITCH OFF STATE 

OUT 

SWON DB(S21) GR1 ! INSERTION LOSS IN dB 

SWOFF DB(S21) GR1A ! LOAD ISOLATION IN dB 

FREQUENCY 

SWEEP 8 22 2 ! FREQUENCY SWEEP FROM 8 TO 22 GHz 

GRID 

GR1 -2 0 .25 

GR1A -20 -80 

OPT 

RANGE 10 20 

LINDON DB(S21 ) = 0 ! MINIMIZE INSERTION LOSS 
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APPENDIX C 

TOUCHSTONE CIRCUIT FILE FOR NODAL VOLTAGE EVALUATION 



This appendix provides the listing of two circuit files using in the Touchstone 
CAD program for the evaluation of nodal voltages by using voltage-controlled voltage 
source circuit model for single pole single throw two-diode switches and single pole 
double throw four-diode switches. 



Cl. SPSTV.CKT single pole single throw nodal voltage evaluation 
C2. SPDTV.CKT single pole double throw nodal voltage evaluation 



Cl. SPSTD.CKT 



! ! SINGLE POLE SINGLE THROW DIODE SWITCH 

! ! NODAL VOLTAGE CALCULATION USING VOLTAGE -CONTROLED VOLTAGE SOURCE 
! ! INPUT PARAMETERS : 

! ! 1. JUNCTION CAPACITANCE, CJ 

! ! 2. JUNCTION RESISTANCE, RJ 

! ! 3. CENTER FREOUENCY, FO 

! ! 4. TERMINAL LOAD IMPEDANCE, RR 

!! 5. TRANSMISSION LINE IMPEDANCE, ZA, ZB 

! ! 6. TRANSMISSION LINE ELECTRICAL LENGTH, EA, EB 

! ! OUTPUT : 

! ! NORMALIZED NODAL VOLTAGE OF THE INTERESTED NODE 
VAR 

! INPUT PARAMETERS: 

FO = 15 ! FREQUENCY IN GHz 

CJ = .15 ! JUNCTION CAPACITANCE IN pF 

RJ = 1.5 ! JUNCTION RESISTANCE IN Ohm 

RR = 50 ! TERMINAL LOAD IMPEDANCE IN Ohm 

! TRANSMISSION LINE IMPEDANCE (Ohm) 

ZA = 36.66 
ZB = 37.30 

! TRANSMISSION LINE ELECTRICAL LENGTH (Degree) 

EA = 108.3 
EB = 58.7 
CKT 

TLIN 1 2 Z A ZA E A EA F A F0 

DEF2P 1 2 SEG1 

TLIN 1 2 Z A ZB E A EB F A F0 

DEF2P 1 2 SEG2 

CAP 1 2 C A CJ 

DEF2P DCJ 

RES 1 2 R A RJ 

DEF2? 1 2 DRJ 

! ON STATE CALCULATION 

5EG1 1 2 

DCJ 2 5 

DRJ 5 0 

SEG2 2 3 

DCJ 3 6 

DRJ 6 0 

SEG 3 4 

RES 4 0 R A RR 

VCVS 4700 M=- . 5 A=0 R1=1E8 R2=0 F=0 T=0 

VCVS 5800 M=-.5 A=0 R1=1E3 R2=0 F=0 T=0 

VCVS 6900 M=-.5 A=0 R1=1E8 R2=0 F=0 T=0 

DEF4P 1789 OUTV1 
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! OFF STATE CALCULATION 
SEG1 11 12 
DRJ 12 0 
SEG2 12 13 
DRJ 13 0 
SEG1 13 14 
RES 14 0 R A RR 
VCVS 14 17 0 0 

VCVS 12 18 0 0 

VCVS 13 19 0 0 

DEF4P 11 17 18 19 OUTV2 
OUT 

! NORMALIZED NODAL VOLTAGE 
OUT VI MAG ( S 2 1 
OUTV1 MAG(S31 
OUT VI MAG ( S41 
OUTV2 MAG (S 21 
OUTV2 MAG(S31 
OUTV2 MAC-(S41 
FREQUENCY 
SWEEP 8 22 2 
GRID 



M=- . 5 A=0 R1=1E8 R2=0 F=0 T=0 
M— .5 A=0 R1 = 1E8 R2=0 F=0 T=0 
M=- . 5 A=0 R1=1E8 R2=0 F=0 T=0 



S21) 


GR1 


S31) 


GR1A 


S41) 


GR1A 


S21 ) 


GR2 


S31) 


GR1A 


S41 ) 


GR2 



! FREQUENCY SWEEP FROM 8 TO 22 GHz 



RANGE 


8 


22 


2 


GR1 


0 


.5 


1 


GR1A 


0 


.15 


.3 


GR2 


0 


.001 


.0001 



C2. SPDTV.CKT 



! ! SINGLE POLE DOUBLE THROW DIODE SWITCH 

! ! NODAL VOLTAGE CALCULATION USING VOLTAGE -CONTROLED VOLTAGE SOURCE 
! ! INPUT PARAMETERS : 

! ! 1. JUNCTION CAPACITANCE, CJ 

! ! 2. JUNCTION RESISTANCE, RJ 

!! 3. CENTER FREQUENCY, FO 

!! 4. TERMINAL LOAD IMPEDANCE, RR 

! ! 5. TRANSMISSION LINE IMPEDANCE, ZA, ZB 

! ! 6. TRANSMISSION LINE ELECTRICAL LENGTH, EA, EB 

! ! OUTPUT : 

! ! NORMALIZED NODAL VOLTAGE OF THE INTERESTED NODE 
VAR 

l INPUT PARAMETERS: 

FO = 15 ! FREQUENCY IN GHz 

CJ = .20 ! JUNCTION CAPACITANCE IN pF 

RJ = 1.5 ! JUNCTION RESISTANCE IN Ohm 

RR = 50 ! TERMINAL LOAD IMPEDANCE IN Ohm 

! TRANSMISSION LINE IMPEDANCE (Ohm) 

ZA = 35.90 
ZB = 47.36 

! TRANSMISSION LINE ELECTRICAL LENGTH ( Degree ) 

EA = 101.20 
Z3 = 34.7 
CKT 

TLIN 1 2 Z A ZA E A EA F A F0 

DEF2P 1 2 SEG1 

TLIN 1 2 Z A ZB E A EB F A F0 

DEF2P 1 2 SEG2 

CAP 1 2 C A CJ 

DEF2P DCJ 

RES 1 2 R A RJ 

DEF2P 1 2 DRJ 

! ON STATE CALCULATION 

SEG1 1 2 
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DCJ 2 3 

DRJ 3 0 

SEG2 2 3 

DCJ 3 9 

DRJ 9 0 

SEG 3 4 

RES 4 0 R A RR 

! OFF STATE CALCULATION 

SEG1 _1 5 

DRJ 3 0 

SEG2 5 6 

DRJ 6 0 

SEG1 6 7 

RES 7 0 R A RR 



vcvs 


4 


10 


0 


0 


M=- . 5 


A=0 


R1=1E8 


R2=0 


F=0 


T=0 


vcvs 


3 


11 


0 


0 


M=- . 5 


A=0 


R1=1E3 


R2=0 


F=0 


T=0 


vcvs 


5 


12 


0 


0 


M=- . 5 


A=0 


R1=1E3 


R2=0 


F=0 


1=0 



DEF4P 1 10 11 12 OUTV1 
OUT 

! NORMALIZED NODAL VOLTAGE 
OUTV1 MAG(S21 ) GR1 
OUT VI MAG(S31) GR1A 
OUTV1 MAG(S41) GR1A 
FREQUENCY 

SWEEP 8 22 2 ! FREQUENCY SWEEP FROM 3 TO 22 GHz 

GRID 

RANGE 8 22 2 

GR1 0.51 
GR1A 0 .5 .1 



80 



APPENDIX D 

FORTRAN PROGRAM FOR DIODE POWER AND HEATING 

CALCULATIONS 

This appendix provides the listing of two Fortran programs written for the 
evaluation of power dissipation and junction temperature in the individual diode of the 
single pole single throw and single pole double diode switches. 

Dl. PSPST diode power dissipation and junction temperature of a single pole 
single throw diode switch. 

D2. PSPDT diode power dissipation and junction temperature of a single pole 
double throw diode switch. 



Dl. PSPST 



5 



10 



15 



20 



CHARACTERS ANSWER 
CHARACTER* 13 FILE 
CHARACTERS SELECT 
CHARACTERS TYPE 
INTEGER I ,NP0RT 
REAL TA,TR,HC,TJON(10) ,TJOFF(10) ,DT 
REALS FRE,Q,CJ,EL(10) , 



RR , PI , PIN , PON , POFF , PLON , PLOFF 



REALS PDON(IO) ,PDOFF(10) \ AMP , VDONM , PEX 

C0MPLEXS6 VIN , ZOUT , Z ( 10 ) , ZINON , ZINOFF , ZIN , I IN , ION , IOFF , ILON , VLON 
COMPLEXS6 ILOFF , VLOFF , VON , VOFF , VS , ZS 

C0MPLEXS6 GPI ( 2 , 2 , 20 ) ,GP (2 , 2 , 20) , GNI (2,2/20) ,GN(2,2,20) 

COMPLEXS6 IDON(lO) .IDOFF(lO) ,VDON(10) 

C0MPLEXS6 VDOFF(IO) 

OPEN (UNIT=2 , FILE= ' OUTPUT 1 ) 

PI = 3.1415926536 

WRITE (6,*) 'POWER DISSIPATION AND JUNCTION TEMPERATURE EVALUATION' 
WRITE (6,*) 'OF SPST DIODE SWITCH ' 

WRITE(6,*) '"PRESS ANY KEY AND RETURN" TO CONTINUE, "Q AND RETURN"TO 
* QUIT' 

READ (5, 10) ANSWER 
FORMAT (Al) 

IF (ANSWER . EQ . ' Q ' ) STOP 

write (6, *) 'Enter' 1 0"to read input from file,"i m to enter inputs inte 
*RACTIVELY’ 

READ (5 , *)NUM 

IF(NUM.NE.O.AND.NUM.NE.l) GOTO 5 

IF (NUM.EQ.O)THEN 

MN = 1 

ELSE 

MN = 5 

END IF 



I F ( NUM . EQ . 0 ) THEN 

WRITES, S 'ENTER THE NAME OF INPUT FILE' 

READ (5 , 15) FILE 

OPEN (UNITS , FILE=FILE ) 

FORMAT(A13) 

END IF 

IF (NUM.EQ . 1 ) THEN 

WRITE (6,*) 'WHAT TYPE OF INPUT SIGNAL "PULSE" OR "CW" ? ' 
WRITE (6,*) 'FOR "PULSE" INPUT ENTER "P" ' 

WRITES,*) 'FOR "CW" INPUT ENTER "C" ' 

WRITE (6,*) 'DO NOT EVALUATE JUNCTION TEMPERATURE ENTER "N" ' 
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25 



30 



35 



40 



END IF 

READ (MN, 11) TYPE 
FORMAT (Al) 

IF (TYPE. NE. 'P 1 .AND. TYPE. NE. 'C' . AND . TYPE . NE . 1 N 1 )GOTO 20 
IF (TYPE. EQ. • N ' )G0T0 30 
IF (TYPE .EQ. 1 C 1 )THEN 

IF (NUM.EQ . 1 ) WRITE (6 , *) 1 ENTER AMBIENT TEMPERATURE (DEGREE C) 1 
READ(MN,*)TA 

IF(NUM.EO.l) WRITE (6,*) 'TOTAL THERMAL RESISTANCE (DEGREE C/WATT) 1 

READ(MN, X )TR 

ELSE 

IF (NUM.EQ • 1)WRITE (6 , *) 1 ENTER AMBIENT TEMPERATURE (DEGREE C) 1 
READ(MN, X )TA 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER THERMAL CAPACITY OF OF I-REGION (MICR 
*0- JOULE/ DEGREE C) 1 
READ(MN,*)HC 

IF (NUM.EQ . 1 ) WRITE (6 , x ) 1 ENTER THE PULSE DURATION (MICRO-SECOND)' 

READ(MN, x )DT 

END IF 

IF (NUM.EQ. 1)THEN 

WRITE (6,*) 1 IS THIS EVALUATION USE THE FIGURE OF MERIT, DIODE Q 1 
WRITE (6,*) 'ANSWER "YES" OR "NO 1 " 

END IF 

READ (MN, 35) SELECT 
FORMAT (A3) 

IF ( SELECT. NE. 'YES' . AND . SELECT .NE . 'NO' ) GOTO 30 
IF ( SELECT . EQ . ' YES ' )THEN 

I F (NUM . EQ . 1 ) WRITE ( 6 , * ) 1 ENTER THE VALUE OF DIODE Q ' 

READ(MN, X )Q 

IF (NUM.EQ. 1) WRITE (6 , *) 1 ENTER THE VALUE OF JUNCTION CAPACITANCE (PF) 



READ (MN , x ) CJ 

IF (NUM.EO . 1 )WRITE (6 , *) 1 ENTER THE EVALUATE FREQUENCY (GHZ)' 

READ(MN,*)FRE 

RR = 1000/ (2*PI*FRE*Q*CJ) 

ELSE 

IF (NUM.EQ. 1)WRITE(6 , *) 1 ENTER THE VALUE OF JUNCTION CAPACITANCE (PF) 

I 



READ(MN,*)CJ 

IF (NUM.EQ . 1)WRITE (6 , *) 1 ENTER THE EVALUATE FREQUENCY (GHZ)' 
READ(MN, X )FRE 

IF(NUM.EQ.1)WRITE(6 ,*) 1 ENTER THE VALUE OF JUNCTION RESISTANCE (OHM 
*) 1 

READ(MN,*)RR 

Q = 1000/(2*PI*FRE*RR*CJ) 

END IF 

VS = (1,0) 

ZS = (50,0) 

IF (NUM.EO . 1 )WRITE (6 , *) 1 ENTER NUMBER OF DIODES IN EACH BRANCH 1 
READ(MN, X )ND 

IF (NUM.EQ . 1 )WRITE (6 , *) 1 ENTER THE EXPECT POWER AT LOAD (WATT) 1 
READ (MN , x )PEX 

IF (NUM . EO . 1 )WRITE (6 , *) 1 ENTER THE OUTPUT IMPEDANCE (COMPLEX) 1 
READ (MN , x ) ZOUT 
NSEG = ND+1 
DO 40 I =1 ,NSEG 

IF (NUM.EO . 1 ) WRITE (6 , x ) 1 TRANSMISSION SEGMENT NUMBER 1 , 1 

IF (NUM.EO . 1 ) WRITE (6 , x ) ' ENTER IMPEDANCE OF TRANSMISSION LINE (COMPL 

*SX) 1 

READ ( MH , * ) Z ( I ) 

I F ( NUM . EQ . 1 ) WRITE ( 6 , * ) 1 ENTER ELECTRICAL LENGTH (DEGREE)' 

RE AD (MN , x ) EL ( I ) 

EL(I) = EL(I)*PI/180 

CONTINUE 

NPORT = ND+NSEG 



C ON STATE CALCULATION 

CALL MATRIX (RR , C J , Z , EL , NSEG , ND , FRE , GP ( 1 , 1 , 1 ) ) 
DO 45 K =1, NPORT 
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IF(K.EQ. 1)THEN 

GPI (1,1,1) = GP (1,1,1) 

GPI (1,2,1) = GP ( 1 , 2 , 1 ) 

GPI (2 ,1,1) = GP ( 2 , 1 , 1 ) 

GPI (2 , 2 , 1 ) = GP ( 2 , 2 , 1 ) 

ELSE 

CALL MATMUL(GPI (1,1 ,K-1 ) ,GP ( 1 , 1 , K) , 2 , 2 , 2 , GPI (1 , 1 ,K) ) 

END IF 

45 CONTINUE 

ZINON = ( GPI ( 1 , 1 , NPORT ) *ZOUT+GPI (1,2, NPORT ) ) / ( GPI ( 2 , 1 , NPORT ) *ZOUT 
#+GPI ( 2 , 2 ,NPORT) ) 



C OFF STATE CALCULATION 

CALL MATRIX(RR, 0 .DO ,Z ,EL ,NSEG,ND , FRE , GN(1 , 1 , 1 ) ) 

DO 50 N =1, NPORT 

IF(N.EQ.1)THEN 

GNI( 1,1,1) = GN ( 1 , 1 , 1 ) 

GNI (1 ,2,1) = GN ( 1 , 2 , 1 ) 

GNI (2,1,1) = GN ( 2 , 1 , 1 ) 

GNI (2 , 2,1) = GN ( 2 , 2 , 1 ) 

ELSE 

CALL MATMUL (GNI (1,1,N-1),GN(1,1,N),2,2, 2 ,GNI ( 1 , 1 ,N) ) 

END IF 

50 CONTINUE 

ZINOFF = ( GNI ( 1 , 1 , NPORT ) *ZOUT+GNI ( 1 ; 2 , NPORT ) ) / (GNI ( 2 , 1 , NPORT ) *ZOUT 
#+GNI (2,2, NPORT ) ) 



C ON-OFF STATE 

VON = ZINON*VS/ (ZINON+VS) 

ION = VON/ZINON 
VOFF = ZINOFF*VS/ (ZINOFF+VS) 

IOFF = VOFF/ZINOFF 
PON = 0.5*REAL(VON*DCONJG(ION)) 

POFF = 0.5*REAL(VOFF*DCONJG(IOFF) ) 

I LON = GPI (1,1, NPORT )*ION-GPI( 2,1, NPORT )*VON 
VLON = GPI (2,2 . NPORT )*VON-GPI ( 1 , 2 ,NPORT)*ION 
PLON = PON-O. 5*REAL(VON*DCONJG( ION) -VLON*DCONJG( ILON) ) 

AMP = PEX/PLON 
WRITE (2, 55) FRE 
WRITE (2, 60) Q 
WRITE (2, 65) ND 
WRITE (2, 70) CJ 
WRITE (2, 75) RR 
WRITE (2, 85) PON*AMP 
WRITE (2,90) POFF*AMP 
WRITE (2, 95) PLON*AMP 

ILOFF = GNI ( 1, 1, NPORT )*IOFF-GNI( 2,1, NPORT )*VOFF 
VLOFF = GNI (2,2 , NPORT )*VOFF-GNI (1,2 , NPORT )*IOFF 
PLOFF = poff-o.5*real(voff*dconjg(ioff)-vloff*dconjg(iloff)) 

WRITE (2, 100) PLOFF*AMP 
55 FORMAT ' 

60 FORMAT 

65 FORMAT 

70 FORMAT 

7 5 FORMAT 

35 FORMAT 

90 FORMAT 

95 FORMAT 

*' ) 

100 FORMAT ( 1 POWER DELIVER TO THE LOAD IN OFF STATE 1 , F10 . 6 , 3X, ' WATT 
*S’ ) 



' , ' OPERATING FREQUENCY 1 , 3X , F5 . 2 , 2X , 1 GHZ 1 ) 

',' FIGURE OF MERIT (DIODE Q)',3X,F8.3) 

', 'TOTAL NUMBER OF DIODE IN EACH BRANCH' , 2X, 12) 

DIODE JUNCTION CAPACITANCE 3X , F5 . 3 , 2X , ' ?F 1 ) 

', 'DIODE JUNCTION RESISTANCE 3X , F5 . 3 , 2X , 1 OHM ' ) 

',' POWER IN ON STATE ' , 2X , F10 . 2 , 3X, 1 WATTS ' ) 

POWER IN OFF STATE' , IX, F10. 2, 3X, 'WATTS' ) 

' 'POWER DELIVER TO THE LOAD IN ON STATE 1 , F10 . 2 , 3X ,' WATTS 



C POWER DISSIPATED IN ON STATE BRANCH 

CALL POWER ( GP ( 1 , 1 , 1 ) , I ON , VON , NPORT , PDON , IDON , VDON ) 
WRITE (2,*) 'POWER DISSIPATION IN DIODE ON STATE' 

DO 105 I = 2 , 2*ND , 2 
WRITE (2 , 115)1/2 ,PDON(I )*AMP 
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105 CONTINUE 

C POWER DISSIPATED IN OFF STATE BRANCH 

CALL POWER (GN ( 1 , 1 , 1 ) , IOFF , VOFF , NPORT , PDOFF , IDOFF , VDOFF ) 

WRITE (2,*) 'POWER DISSIPATION IN DIODE OFF STATE' 

DO 110 J = 2 , 2*ND , 2 
WRITE (2 , 115) J/2 , PDOFF (J)*AMP 
110 CONTINUE 

115 FORMAT (' ', 'POWER DISSIPATED IN DIODE NO . ' , IX , 12 , 2X, F10 . 5 , 2X , ' WATT 

*S' ) 

C JUNCTION TEMPERATURE CALCULATION 

C JUNCTION TEMPERATURE IN ON STATE BRANCH 

IF (TYPE .EQ . ' N ' )STOP 
IF(TYPE.EQ. 'C' )THEN 
CALL THERM 1 ( ND , T JON , PDON , AMP , TA , TR ) 

ELSE 

CALL THERM2 ( ND , T JON , PDON , AMP , TA , DT , HC ) 

END IF 

WRITE (2,*) 'JUNCTION TEMPERATURE ON STATE BRANCH' 

DO 120 I = 1 ,ND 
WRITE (2 ,82)1 ,TJON(I ) 

120 CONTINUE 

C JUNCTION TEMPERATURE IN OFF STATE BRANCH 

WRITE (2,*) 'JUNCTION TEMPERATURE OFF STATE BRANCH' 

IF (TYPE. EO. ' C' )THEN 

CALL THERM1 (ND , TJOFF , PDOFF , AMP , TA , TR ) 

ELSE 

CALL THERM2 (ND , TJOFF , PDOFF , AMP , TA , DT , HC ) 

END IF 

DO 125 J = 1 , ND 
WRITE (2, 130 )J, TJOFF (J) 

125 CONTINUE 

130 FORMAT ( ' ',' JUNCTION TEMPERATURE DIODE NO . ' , 2X, 12 , 2X, F6 . 2 , 2X, ' DEGR 

*EE C ' ) 

STOP 

END 



c SUBROUTINE MATRIX 

SUBROUTINE MATRIX ( RR , C J , Z , EL , NSEG , ND , FRE , GP ) 

REAL *8 RR , CJ , EL ( 10 ) , FRE , PI 
C0MPLEX*16 2(10), 21 

COMPLEX*16 G(2,2,20) ,G1(2,2,20) ,GP(2,2,20) 

PI = 3.1415926536 
DO 10 K = 1 ,NSEG 

G( 1 , 1 , K) = DCMPLX(COS (EL(K) ) .0 .DO) 

G(1,2,K) = DCMPLX(0.D0,1.D0)*z(k)*DCMPLX(SIN(EL(K)),0.D0) 
G(2,1,K) = DCMPLX(0.D0,1.D0)/Z(K)*DCMPLX(SIN(EL(K) ) ,0.D0) 
G(2 , 2 ,K) = G(1 , 1 ,K) 

10 CONTINUE 

DO 20 L=1 ,ND 
IF(CJ.EQ.0.D0)THEN 
Z1 = DCMPLX ( RR , 0 . DO ) 

ELSE 

Z1 = DCMPLX (RR , -1000/ (2*PI*FRE*CJ) ) 

END IF 

C-i( 1,1,1) = DCMPLX(1 .DO ,0 .DO) 

G1(1,2,L) = DCMPLX (0 .DO , 0 . DO ) 

G1 (2 , 1 ,L)= 1/21 

G1 ( 2 , 2 , L) = G1 (1 , 1 , L) 

20 CONTINUE 

DO 30 M =1 ,NSEG 
N = 2*M-1 

GP ( 1 , 1 , N ) = G( 1 , 1 ,M) 

GP( 1 , 2 ,N) = G(1,2,M) 

GP ( 2 , 1 ,N) = G(2,1,M) 

GP ( 2 , 2 , N ) = G(2,2,M) 

30 CONTINUE 

DO 40 I =1 ,ND 
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40 



J = 2*1 

GP ( 1 , 1 , J ) = G1 (1 , 1 , 1 ) 

GP(1,2,J) = G1 ( 1 , 2 , 1 ) 

GP ( 2 , 1 , J ) = G1 ( 2 , 1 , 1 ) 

GP ( 2 , 2 , J ) = G1 (2 , 2 , 1 ) 

CONTINUE 

RETURN 

END 



C SUBROUTINE MATMUL 

SUBROUTINE MATMUL ( AA , BB , NN , MM , RR , CC ) 

INTEGER NN , MM , RR 

COMPLEX*16 AA(2,2,1),BB(2,2,1),CC(2,2,1) 

DO 10 I = 1 ,NN 
DO 20 J = 1,RR 
CC(I , J, 1) = 0 
DO 30 K =1 ,MM 

CC(I,J,1) = CC(I,J,1) + AA ( I , K , 1 ) *BB ( K , J , 1 ) 
30 CONTINUE 

20 CONTINUE 

10 CONTINUE 

RETURN 
END 



c SUBROUTINE POWER 

SUBROUTINE POWER ( G , I , V , N , PD , I ID , VID ) 

REAL*8 PD (10) 

COMPLEX*16 G(2,2,20) ,V,I,IID(20) , VID (20) , I0D(20) ,VOD(20) 

DO 10 J=2 ,N , 2 

IID(J) = G(1 ,1 , J-1)*I - G(2, 1 , J-1)*V 
VID(J) = G( 2 , 2 , J-l )*V - G(1 , 2 , J-1)*I 
IOD(J) = G(1,1,J)*IID(J) - G(2 , 1 , J)*VID( J) 

VOD(J) = G(2 .2 , J)*VID(J) - G( 1 , 2 , J)*IID( J) 

PD ( J) = 0.5*REAL(VID(J)*DCONJG(IID(J))-VOD(J)*DCONJG(IOD(J))) 
I = IOD(J) 

V = VOD(J) 

10 CONTINUE 
RETURN 
END 



c SUBROUTINE THERM1 

SUBROUTINE THERM1 (ND , T J , PD , AMP , TA , TR ) 
INTEGER ND 

REAL TJ(10) ,AMP,TA,TR 
REAL*8 PD (10) 

DO 5 1=1, ND 

TJ(I ) = TA + ( PD ( 1*2 ) *AMP ) *TR 
5 CONTINUE 

RETURN 
END 



c SUBROUTINE THERM2 

SUBROUTINE THERM2 (ND , T J , ?D , AMP , TA , DT , HC) 
INTEGER ND 

REAL T J ( 1 0 ) , AMP , TA , DT , HC 
REAL*8 PD (10) 

DO 5 J =1,ND 

TJ( J) = TA + (PD(J*2)*AMP)*DT/HC 
5 CONTINUE 

RETURN 
END 
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D2.PSPDT 



CHARACTER*! ANSWER 
CHARACTER* 13 FILE 
CHARACTER*3 SELECT 
CHARACTER*! TYPE 
INTEGER I , NPORT 

REAL TA,TR,HC,TJON(10) ,TJOFF(10) ,DT 

REAL*3 FRE,Q,CJ,EL(10) ,RR,PI ,PIN,PON,POFF,PLON,PLOFF 
REAL *8 PDON ( 10 ) , PDOFF ( 10 ) , AMP , VDONM , PEX 

COMPLEX* 16 VIN , ZOUT , Z ( 10 ) , ZINON , ZINOFF , ZIN , I IN , ION , IOFF , ILON , VLON 
COMPLEX*16 ILOFF , VLOFF , VON , VOFF , VS , ZS 

COMPLEX*16 GPI (2,2,20) ,GP(2,2,20) ,GNI(2 ,2 ,20) ,GN(2 ,2 ,20) 

COMPLEX*16 IDON(IO) .IDOFF(IO) ,VDON(10) 

COMPLEX*16 VDOFF(IO) 

OPEN(UNIT=2 , FILE= 1 OUTPUT 1 ) 

PI = 3.1415926536 

5 WRITE (6,*) 'POWER DISSIPATION AND JUNCTION TEMPERATURE EVALUATION' 

WRITE (6,*) 'OF SPDT DIODE SWITCH ' 

WRITE (6,*) "'PRESS ANY KEY AND RETURN" TO CONTINUE, "Q AND RETURN"TO 
* QUIT' 

READ(5. 10) ANSWER 
10 FORMAT (Al) 

IF ( ANSWER. EQ. • Q ' )STOP 

WRITE (6 , *) 1 cNTER"0"TO READ INPUT FROM FILE,"l"TO ENTER INPUTS INTE 
*RACTIVELY ' 

READ ( 5 , * ) NUM 

IF (NUM .NE . 0 . AND .NUM.NE . 1 ) GOTO 5 

IF (NUM.EQ . 0)THEN 

MN = 1 

ELSE 

MN = 5 

END IF 

IF (NUM.EQ . 0 ) THEN 

WRITE (6,*) 'ENTER THE NAME OF INPUT FILE' 

READ (5 , 15) FILE 
OPEN(UNIT=l,FILE=FILE) 

15 FORMAT (A13) 

END IF 

IF(NUM.EQ.1)THEN 

20 WRITE (6,*) 'WHAT TYPE OF INPUT SIGNAL "PULSE" OR "CW" ? ' 

WRITE (6, *)' FOR "PULSE" INPUT ENTER "P" 1 
WRITE (6,*) 'FOR "CW" INPUT ENTER "C" ' 

WRITE (6 , *) ' DO NOT EVALUATE JUNCTION TEMPERATURE ENTER "N" ' 

END IF 

READ (MN, 11) TYPE 
25 FORMAT (Al) 

IF (TYPE .NE . ' P ' .AND. TYPE. NE. 'C' .AND .TYPE .NE .' N ' )GOTO 20 
IF (TYPE. EO. ' N ' )GOTO 30 
I F ( TYPE . EO . ' C 1 ) THEM 

IF (NUM. EC 7l ) WRITE (6 , *) ' ENTER AMBIENT TEMPERATURE (DEGREE C) ' 

READ ( MN , * ) TA 

IF (NUM. EO . 1 ) WRITE (6 , *) 1 TOTAL THERMAL RESISTANCE (DEGREE C/WATT) ' 

READ(MN,*)TR 

ELSE 

IF (NUM.EQ . 1 ) WRITE (6 , *) ' ENTER AMBIENT TEMPERATURE (DEGREE C) ' 

read(mn,*}ta 

IF(NUM .EQ . 1 )WRITE (6 , *) 'ENTER THERMAL CAPACITY OF OF I-REGION (MICR 
*0- JOULE/ DEGREE C) ' 

READ(MN,*)HC 

IF(NUM.EQ . 1 )WRITE (6 , *) 'ENTER THE PULSE DURATION (MICRO-SECOND)' 

READ(MN,*)DT 

END IF 
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30 



35 



40 



I F (NUM . EQ . 1 ) THEN 

WRITE(6 ,*) ' IS THIS EVALUATION USE THE FIGURE OF MERIT , DIODE O' 

WRITE (6,*) 'ANSWER "YES" OR "NO" 1 
END IF 

READ(MN, 35) SELECT 
FORMAT (A3) 

IF(SELECT.NE. ‘YES' . AND . SELECT .NE . 1 NO 1 )GOTO 30 
IF ( SELECT . EO . 'YES 1 )THEN 

IF(NUM.E0.1)WRITE(6,*) 'ENTER THE VALUE OF DIODE 0 1 
READ(MN,*)Q 

IF(NUM. EQ . 1 ) WRITE (6 , *) ' ENTER THE VALUE OF JUNCTION CAPACITANCE (PF) 
READ(MN,*)CJ 

* i 

IF (NUM.EQ . 1 )WRITE (6 , *) ' ENTER THE EVALUATE FREQUENCY (GHZ)' 

READ(MN,*)FRE 

RR = 1000/(2*?I*FRE*0*CJ) 

ELSE ~ 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER THE VALUE OF JUNCTION CAPACITANCE (PF) 

* i 



READ(MN,*)CJ 

IF(NUM'.EQ.1)WRITE(6,*) 'ENTER THE EVALUATE FREQUENCY (GHZ)' 
READ(MN,*)FRE 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER THE VALUE OF JUNCTION RESISTANCE (OHM 
') 1 

READ (MN , *)RR 

Q = 1000/ (2*?I*FRE*RR*CJ ) 

END IF 
VS = (1,0) 

ZS = (50,0) 

IF (NUM.EQ . 1 ) WRITE (6 , *) 1 ENTER NUMBER OF DIODES IN EACH BRANCH 1 
READ (MN , x ) ND 

IF (NUM. EQ.l) WRITE (6,*) 'ENTER THE EXPECT POWER AT LOAD (WATT)' 
READ(MN,*)PEX 

IF (NUM. EQ.l) WRITE (6,*) 'ENTER THE OUTPUT IMPEDANCE (COMPLEX)' 

READ (MN , * j ZOUT 
NSEG = ND+1 
DO 40 I =1,NSEG 

IF(NUM. EQ . 1 ) WRITE (6 , * j ' TRANSMISSION SEGMENT NUMBER', I 



IF (NUM. EQ.l) WRITE (6 
*EX) ' 

READ(MN,*)Z(I) 



'ENTER IMPEDANCE OF TRANSMISSION LINE (COMPL 



IF (NUM.EQ. 1 )WRITE (6 , *) ' ENTER ELECTRICAL LENGTH (DEGREE) 
READ(MN,*)EI 



. . : )EL(I) 

EL(l) = EL(I)*PI/180 

CONTINUE 

NPORT = ND+NSEG 



C 



45 



ON STATE CALCULATION 

CALL MATRIX(RR,CJ,Z,EL, NSEG, ND,FRE,GP(1, 1,1)) 

DO 45 K =1, NPORT 

IF(K.EQ.1)THEN 

GPI( 1,1,1) = GP( 1,1,1) 

GPI (1,2,1) = GP ( 1 , 2 , 1 ) 

GPI (2,1,1) = GP ( 2 , 1 , 1 ) 

GPI (2,2,1) = GP ( 2 , 2 , 1 ) 

ELSE 



CALL MATMUL ( G? I ( 
END IF 



K-l) , 3? ( 1 , 1 , X ) ,2,2,2,G?I(1,1,K)) 



CONTINUE 

ZINON = ( GP I ( 1 , 1 , NPORT ) *ZOUT+GP 1(1,2, NPORT ) ) / ( GPI ( 2 , 1 , NPORT ) *ZOUT 
#+GPI (2, 2, NPORT) ) 



C OFF STATE CALCULATION 

CALL MATRIX (RR, 0 .DO , Z ,EL ,NSEG,ND , FRE ,GN( 1 , 1 , 1 ) ) 
DO 50 N =1, NPORT 
IF (N. EQ.l) THEN 
GNI (l,l,l) = GN ( 1 , 1 , 1 ) 

GNI (l,2,l) = GN ( 1 , 2 , 1 ) 



87 



GNI (2,1,1) = GN(2,1,1) 

GNI (2,2,1) = GN ( 2 , 2 , 1 ) 

ELSE 

CALL MATMUL(GNI(1,1,N-1) ,GN(1,1,N) , 2 , 2 , 2 , GNI ( 1 , 1 ,N) ) 

END IF 

50 CONTINUE 

ZINOFF = ( GNI ( 1 , 1 , NPORT ) *ZOUT+GNI (1,2, NPORT ) ) / (GNI ( 2 , 1 , NPORT ) *ZOUT 
#+GNI (2,2 ,NPORT ) ) 



ON-OFF 
ZIN = 
VIN = 

I IN = 
PIN = 
ION = 
IOFF = 
PON = 
POFF = 
ILON = 
VLON = 
PLON = 
AMP 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE ( 



STATE 

ZINON*ZINOFF/ (ZINON+ZINOFF) 

ZIN*VS/ (ZIN+ZS) 

VIN/ZIN 

REAL( .5*(VIN*DCONJG(IIN))) 

IIN*ZINOFF/ (ZINON+ZINOFF) 

IIN - ION 

0 .5*REAL(VIN*DCONJG(ION)) 

0 . 5*REAL (VIN*DCONJG (IOFF ) ) 

GPI (1,1, NPORT ) *ION-GPI (2,1, NPORT ) *VIN 

GPI (2,2. NPORT ) *VIN-GPI (1,2, NPORT ) *ION 

PON-O . 5^REAL ( VIN*DCONJG (ION) -VLON*DCONJG (ILON) ) 

PEX/PLON 



55 


ILOFF = 
VLOFF = 
PLOFF = 
WRITE (2 
FORMAT ( 


60 


FORMAT ( 


65 


FORMAT ( 


70 


FORMAT ( 


75 


FORMAT ( 


80 


FORMAT ( 


85 


FORMAT ( 


90 


FORMAT ( 


95 


FORMAT ( 


100 


*') 

FORMAT ( 




*S' ) 



2,55 

2,60 

2,55 

2,70 

2,75 

2,30 

2,85 

2,90' 

2,95:' 



FRE 

Q 

ND 

CJ 

RR 

PIN* AMP 
PON*AMP 
?OFF*AMP 
PLON*AMP 



, 'OPERATING FREQUENCY 1 , 3X, F5 . 2 , 2X, ' GHZ ' ) 

, 1 FIGURE OF MERIT (DIODE Q)',3X,F8.3) 

, 'TOTAL NUMBER OF DIODE IN EACH BRANCH 1 ,2X, 12) 

, 1 DIODE JUNCTION CAPACITANCE 1 , 3X, F5 . 3 , 2X , 1 PF ' ) 

, 1 DIODE JUNCTION RESISTANCE ', 3X.F5 . 3 , 2X, ' OHM ' ) 

, 1 INPUT POWER' ,F10. 2, 3X, 'WATTS' ) 

,' POWER IN ON STATE BRANCH' ,2X,F10 .2, 3X, 'WATTS ' ) 

, 'POWER IN OFF STATE BRANCH ', IX, F10 . 2 , 3X, ' WATTS 1 ) 

, 'POWER DELIVER TO THE LOAD IN ON STATE ', F10 . 2 , 3X, ' WATTS 

, 'POWER DELIVER TO THE LOAD IN OFF STATE ', F10 .6 , 3X, ' WATT 



C POWER DISSIPATED IN ON STATE BRANCH 

VON = VIN 

CALL POWER ( GP ( 1 , 1 , 1 ) , ION , VON , NPORT , PDON , IDON , VDON) 

WRITE (2,*) 'POWER DISSIPATION IN DIODE ON STATE' 

DO 105 I = 2 , 2*ND , 2 
WRITE ( 2 , 1 1 5 ) I / 2 , PDON ( I ) *AMP 
105 CONTINUE 

C POWER DISSIPATED IN OFF STATE BRANCH 

VOFF = VIN 

CALL POWER (GN ( 1 , 1 , 1 ) , IOFF , VOFF , NPORT , PDOFF , IDOFF , VDOFF ) 

WRITE (2,*) 'POWER DISSIPATION IN DIODE OFF STATE' 

DO 110 J = 2 , 2*ND , 2 
WRITE (2, 115 )J/2, PDOFF (J)*AMP 
110 CONTINUE 

115 FORMAT ( ' 1 , 'POWER DISSIPATED IN DIODE NO . ' , IX, 12 , 2X, F10 . 5 , 2X, ' WATT 
*S') 

C JUNCTION TEMPERATURE CALCULATION 

C JUNCTION TEMPERATURE IN ON STATE BRANCH 
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IF (TYPE. EQ. ' N 1 )STOP 
IF(TYPE.EQ. 1 C 1 ) THEN 
CALL THERM1 (ND , T JON , PDON , AMP , TA , TR ) 

ELSE 

CALL THERM2 (ND , T JON , PDON , AMP , TA , DT , HC ) 

END IF 

WRITE (2,*) 'JUNCTION TEMPERATURE ON STATE BRANCH 1 
DO 120 I = 1 , ND 
WRITE (2,82)1, T JON ( I ) 

120 CONTINUE 

C JUNCTION TEMPERATURE IN OFF STATE BRANCH 

WRITE (2 , *) ' JUNCTION TEMPERATURE OFF STATE BRANCH' 

IF (TYPE. EQ. ' C ' )THEN’ 

CALL THERM1 (ND , TJOFF , PDOFF , AMP , TA , TR) 

ELSE 

CALL THERM 2 (ND , TJOFF , PDOFF , AMP , TA , DT , HC ) 

END IF 

DO 125 J = 1 ,ND 
WRITE(2 , 130) J,TJOFF(J) 

125 CONTINUE 

130 FORMAT ( ' JUNCTION TEMPERATURE DIODE NO . ' , 2K, 12 , 2X, F6 . 2 , 2X, 1 DEGR 

*EE C') 

STOP 

END 



c 



10 



20 



30 



40 



SUBROUTINE MATRIX 

SUBROUTINE MATRIX (RR , C J , Z , EL , NSEG , ND , FRE , GP ) 

REAL*8 RR,CJ,EL(10) , FRE, PI 
COMPLEX*16 Z(10),Z1 

COMPLEX*16 G(2,2,20) ,G1(2,2,20) ,GP(2,2,20) 

PI = 3.1415926536 
DO 10 K = 1 ,NSEG 

G(1 , 1 ,K) = DCMPLX(COS (EL(K)) . 0 .DO ) 

G( 1 , 2 ,K) = DCMPLX(0.D0,1.D0)*Z(K)*DCMPLX(SIN(EL(K)) 
G(2 , 1 ,K) = DCMPLX(0.D0,1.D0)/Z(K)*DCMPLX(SIN(EL(k)) 
G(2 , 2 ,K) = G(1,1,K) 

CONTINUE 
DO 20 L=1,ND 
IF(CJ.EQ.O .DO) THEN 
Z1 = DCMPLX (RR , 0 . DO ) 

ELSE 

Z1 = DCMPLX (RR , -1000/ (2*PI*FRE*CJ ) ) 

END IF 

G1 ( 1 , 1 , L) = DCMPLX ( 1 . DO , 0 . DO ) 

G1(1,2,L) = DCMPLX ( 0 . DO , 0 . DO ) 

G1 (2, 1 ,L)= 1/Z1 

G1 (2 , 2 ,L) = G1 (1 , 1 ,L) 

CONTINUE 

DO 30 M =1 ,NSEG 

N = 2*M-1 



GP ( 


1,1, N) 


= G( 


1,1, M 


GP < 


1 , 2 ,N) 


= G < 


1,2,M 


GP! 


2 , 1 ,N) 


= G ( 


2 , 1 ,M) 


GP! 


2 , 2 ,N ) 


= Gi 


2 , 2 , M ) 



CONTINUE 



DO 40 I = 1 , ND 
J - 2^1 



OP ( 1 , 1 , J 
GP(1 , 2 , J 
GP ( 2 , 1 , J 
GP(2 , 2 , J 
CONTINUE 
RETURN 
END 



) = Gl( 




) = Gil 


1,2,1 


) = Gil 


2,1,1 


) = Gil 


;2,2,i) 



0.D0) 
0 .DO) 
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c SUBROUTINE MATMUL 

SUBROUTINE MATMUL ( AA , BB , NN , MM , RR , CC ) 

INTEGER NN , MM , RR 

COMPLEX*16 AA(2,2,1) , BB ( 2 , 2 , 1 ) , CC ( 2 , 2 , 1 ) 

DO 10 I = 1 , NN 
DO 20 J = 1,RR 
CC(I,J,1) = 0 
DO 30 K =1 ,MM 

CC(I,J,1) = CC(I,J,1) + AA(I ,K, 1 )*BB (K, J , 1 ) 
30 CONTINUE 

20 CONTINUE 

10 CONTINUE 

RETURN 
END 



c SUBROUTINE POWER 

SUBROUTINE POWER (G , I , V , N , PD , IID , VID ) 

REAL*8 PD (10) 

C0MPLEX*16 G(2,2,20) ,V,I,IID(20) , VID (20) ,IOD(20) ,VOD(20) 

DO 10 J=2 ,N, 2 

IID(J) = G(1,1,J-1)*I - G(2 , 1 , J-l )*V 
VID(J) = G(2,2,J-l)*V - G( 1 , 2 , J-l ) *1 
IOD(J) = G(l,l, J)*IID(J) - G(2,1,J)*VID(J) 

VOD(J) = G(2 .2, J)*VID(J) - G( 1 , 2 , J) *IID ( J) 

PD(J) = 0.5*REAL(VID(J)*DCONJG(IID(J))-VOD(J)*DCONJG(IOD(J))) 
I = IOD(J) 

V = VOD(J) 

10 CONTINUE 
RETURN 
END 



c SUBROUTINE THERM1 

SUBROUTINE THERM1 (ND , TJ , PD , AMP , TA , TR) 
INTEGER ND 

REAL T J ( 1 0 ) , AMP , TA , TR 
REAL*8 PD (10) 

DO 5 1=1, ND 

TJ (I ) = TA + (PD(I*2)*AMP)*TR 
5 CONTINUE 

RETURN 
END 



c SUBROUTINE THERM 2 

SUB ROUT INE THERM2 ( ND , T J , PD , AMP , TA , DT , HC ) 
INTEGER ND 

REAL T J ( 1 0 ) , AMP , TA , DT , HC 
REAL*8 PD (10) 

DO 5 J =1 ,ND 

TJ(J) = TA + (PD(J*2)*AMP)*DT/HC 
5 CONTINUE 
RETURN 
END 
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